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CN ; ABSTRACT 

We present results from the observations and modeling of seventeen Class I cores with the Submil- 
limetre Common Users Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT). Our 
sample consists of cores with 64< T;,o; (K)<270, 0.2 < Lofcs/L© < 12, and 50 < Lobs/Lsmm < 1000. By 
QQ ' modeling the transfer of radiation through the envelope for nine cores, we find, for a power law distri- 

bution n{r) — nf{r/rf)~P, the average and standard deviation (p) = 1.6 ± 0.4 and a median of p = 1.8. 
• However, the inclusion of a disk or other point-like component can cause the derived p to be shallower 

' by as much as 0.5. We discuss uncertainties due to the interstellar radiation field (ISRF), disks, dust 

opacity, and outer radii in our modeling results. We find no evidence for a truncated outer radius or 
radially variant dust properties in most sources. Uncertainty in the strength of the ISRF and possible 
existence of a disk contribute the greatest uncertainty in p. In addition, we test the Shu collapse model 

■ for our sources and discuss the application of simpler analyses that derive a density power law distribu- 
' tion directly from the slope of the intensity radial profile. The total mass of the envelope in our sample 
, has a range of 0.04 < Me„„/M0 < 5.0, but these masses disagree with the virial masses derived from 
' molecular line observations, indicating that observations of molecular lines do not trace the mass in some 
I Class I cores. We also discuss several sources individually. In particular, IRAS 03256-1-3055, with its 

'p^, unique morphology, is an ideal object for testing theories of fragmentation in the formation of low-mass 

Q ' protostars. Also, we note the possibility, through some simple calculations, that IRAS 043854-2550 is 

^ I a young, forming substellar object. Finally, we discuss the nature of these sources in light of various 

C/3 . evolutionary indicators and find that Tboi and Lobs/Lsmm are often inconsistent in distinguishing Class 

' from Class I cores. We note that, in this sample, the Lobs/Lsmm criterion redefines many of these Class 

^ , I sources (by Tboi) as Class sources. 
• 1— ( . 

k> , Subject headings: stars: formation, low-mass — ISM: dust — ISM: individual (CB230, L1251B, IRAS 

■ 04166-1-2706, IRAS 04169-f 2702, IRAS 04239-1-2436, IRAS 04248+2612, IRAS 
C3 ■ 04295-1-2251, IRAS 04361-^2547, IRAS 04381-^2540, IRAS 04385-^2550, IRAS 

03256-1-3055) 
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1. INTRODUCTION 

The formation of low-mass stars has received consider- 
able attention, but the process by which material moves 
from the outer envelope to the central protostar is still de- 
bated. Theories differ, primarily, in their prediction of the 
density, n(r), and velocity, v{r), as a function of radius, 
r. Larson (1969) and Penston (1969) proposed a model 
with initial conditions of constant density and zero veloc- 
ity. Their model encounters problems due to its prediction 
of supersonic infall velocities that are not observed in re- 
gions of low-mass star formation (Zhou 1992). Shu (1977) 
proposed a scenario where infall begins from the inside, 
i.e., inside-out collapse. The initial density configuration 
is a singular isothermal sphere (SIS) with n(r) oc . 
Collapse begins at the center and propagates outward at 
the sound speed, Oe//, leaving behind a density distribu- 
tion that approaches n(r) oc r~^-^ . Some perturbations to 
the Shu model include rotation (Terebey, Shu, & Cassen 



1984) and magnetic fields (GaUi & Shu 1993). McLaugh- 
lin & Pudritz (1997) proposed a different model in which 
accretion increases with time. Their scenario was similar 
to Shu (1977) — a collapse front propagating into a static 
envelope — but, because of differing initial conditions and 
a polytropic equation of state, they predict a much shal- 
lower density profile in the static envelope {n oc r^^). Fos- 
ter & Chevalier (1993) conducted hydrodynamical calcu- 
lations of the collapse of a Bonnor-Ebert sphere and found 
the density approaches a r^^ distribution at small radii. 
Other models exist with various predictions in the density 
distribution, n oc r^^. Therefore, observational determi- 
nations of this density distribution can help us to discern 
the physical processes governing the formation of low- mass 
stars. 

Classification of the stages in protostellar evolution has 
grown more detailed. Initially, Lada et al. (1987) used the 
infrared spectral index to define Classes I-III. Later, Myers 
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& Ladd (1993) proposed a classification of these forming 
stars based on their bolometric temperature (Tboi)', this 
classification scheme was further codified by Chen et al. 
(1995). Andre et al. (1993) defined a criterion to discern 
Class from Class I cores; the ratio of the submillime- 
ter to bolometric luminosity. In this paper, we refer to 
the inverse of this quantity {Lboi/Lsmm), a quantity that 
presumably increases with time. 

The luminosity of a low-mass core (protostar+envelope) 
is emitted mostly at infrared wavelengths due to the re- 
processing of higher energy radiation from the central 
protostar by dust grains in the envelope. These heated 
grains also emit detectable radiation at submillimcter 
wavelengths. Because the cloud is optically thin to this ra- 
diation, a measure of this dust emission renders a sampling 
of the entire mass in the cloud. In the past, astronomers 
have used single bolometers to detect this emission, but 
resolution and sensitivity limited their progress. Submil- 
limcter bolometer arrays, such as SCUBA, have allowed us 
to map this dust emission from the protostellar envelopes 
and, for the first time, draw conclusions about the nature 
and distribution of mass within these envelopes. 

In this paper, we present the observations and modeling 
of Class I cores. We describe our sample selection, ob- 
servations, data reduction, and some basic data analysis 
in Section 2. In Section 3, we discuss various aspects of 
modeling these sources, power law models were the pri- 
mary focus, but wc; also present simpler top hat models 
of the density distribution. Also, we report on Shu col- 
lapse models for nine sources. Additionally, in this third 
section, wc compare the application of simpler analyses to 
our more detailed models. Also, we evaluate some pos- 
sible scenarios where a disk is present and how the disk 
affects the modeled density profile. In Section 4, we dis- 
cuss some individual sources. Section 5 holds a discussion 
of masses, classification, and the density structure of Class 
I protostars. Finally, in Section 6, we summarize our con- 
clusions. The frugal reader, who wishes only for our major 
conclusions, can focus primarily on Sections 3, 5, and 6. 

2. SCUBA OBSERVATIONS 
2.1. Sample 

In Shirley et al. (2000, hereafter, Paper I), we presented 
the results from observations of a sample ranging from 
Class —1 to Class I cores. In that sample, the presence 
of multiple sources within the observed field was common 
for the Class I sources but uncommon for the Class ob- 
jects. Therefore, Shirley et al. suggested that many Class 
I sources might appear as such because of the existence of 
miiltiple sources within the IRAS beam. This suggestion 
was one motivation for observing this sample of Class I 
objects: to determine the existence and effect of multiplic- 
ity. In addition, we hoped to present a similar treatment 
and analysis as in Papers I-III (Shirley et al., 2000; Evans 
et al. 2001; Shirley et al. 2002) —i.e., determining source 
characteristics and proposing possible density and temper- 
ature profiles through simulations with a one-dimensional 
radiative transfer code. 

We have selected a sample of sixteen sources; these 
are listed in Table 1 with their positions in B1950.0 and 
J2000.0 coordinates. The positions were obtained, mostly, 
from the IRAS PSC as well as from 2 /xm observations 



pursued by Tamura et al. (1991) and Moriarty-Schieven 
et al. (1994). 

The distances to these sources range from 140 to 450 
parsecs (pc), but most of the objects (13) are in the nearby 
Taurus star- forming region at 140 pc (Elias 1978). We 
adopt a distance of 450 pc for CB230 (Launhardt & Hen- 
ning 1997), 300 pc for L1251B (Kun & Prusti 1993), and 
320 pc for IRAS 03256-1-3055 in Perseus (de Zeeuw 1999). 
The distance for Perseus is based on recent Hipparcos ob- 
servations and differs from that used in Papers I & III (220 
pc). 

We classify the cores in our sample by their bolometric 

temperature {Tboi > 70 K for Class I). Therefore, two ob- 
jects, L1251B and IRAS 03256-^3055 (with Tboi < 70 K), 
are not inc;luded in the average values for various param- 
eters describing Class I cores. 

2.2. Observations 

We observed this sample of Class I sources with SCUBA, 
the Submillimcter Continuum Users Bolometer Array 
(Holland et al. 1999), on the James Clerk Maxwell Tele- 
scope (JCMT) atop Mauna Kea, Hawaii. Over four nights 
in November 1999 and February 2000, we mapped each 
of these cores at 850 and 450 /im by the 64-point jiggle 
map method that fully samples a field of 2.'3. Since we 
expected to see very extended emission, we adopted a five- 
point mapping technique with 30" spacings between each 
of the five 64-point jiggle maps. By this method, we re- 
moved the effects of bad bolometers and integrated for a 
longer time on the central 2' of the map. In some sources, 
we increased the coverage of our maps to include extended 
emission. 

We performed pointing and skydip measurements be- 
tween each five-point map. The average pointing RMS, 
over the four nights of observations, was ±2'.'0. We mea- 
sured the sky opacity via the skydip method at 450 and 
850 /im ((t45o) = 1.07 and (t85o) = 0.21) while also record- 
ing the opacity as measured by the CSO radiometer at 225 
GHz {{tcso) = 0.05). In our calculation of ts5o and T450, 
we used the updated hot and cold load temperatures as 
given by Archibald et al. (2000). In addition, we find a 
relationship between rgso and tcso that is in agreement 
with that presented by Archibald et al. (2000). We list all 
opacity measurements in Table 2. 

Observations of Uranus (November 1999) and AFGL 
618 (February 2000) were used to determine the beam pro- 
file of the .JCMT. In Figure 1, we show the radial inten- 
sity profiles of these point sources for various times dur- 
ing the night — our observing shift was from about 1730- 
0130, Hawaii- Aleutian Standard Time (HST). Two points 
should be considered in regard to these beam profiles: the 
significant and variable nature of the sidelobes and the 
changing FWHM size of the beam {9mb)- First, Mueller 
et al. (2002) noted the significance of sidelobes in their 
treatment of dust continuum emission from massive star- 
forming regions. They find that disregard of these side- 
lobes and the assumption of a Gaussian beam can cause 
the derived density power law index to be shallower by 
~0.4. Indeed, the beam profiles in Figure 1 show consid- 
erable power (about 1/10 of the maximum) in the second 
and, sometimes, third sidelobes. Additionally, we find the 
sidelobes become smaller as the night progresses and the 
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dish cools. Second, for both nights, the 450 /im beam pro- 
file has a large FWHM (~11") at about 1800 HST, but 
the beamsize decreases to 6mb ~ 8" at about midnight 
(HST). For November 1999, as shown in Figure 1, the 
450 /im FWHM decreased by 3" as the night progressed. 
However, because the S/N is much better and the beam 
more stable at 850 /um, we relied more heavily on it in the 
determination of certain source properties (e.g., the decon- 
volved FWHM size). The 850 /im beam profile is also vari- 
able {SFWHM ~ 1 — 2"), however, and we propagate this 
variability as an uncertainty in the measured properties de- 
pendent on the beamshape. Additionally, in our modeling, 
we used the observed beam profile that is closest in time to 
the source observation. In one case (CB230), where there 
were no beam profiles within one hour in time, we took 
the time-weighted average of two neighboring beam pro- 
files; in all other cases, the beam profiles were taken from 
observations of Uranus or AFGL 618 that were within one- 
half hour of the source observation. While the existence 
of and variations in the sidelobes significantly affect the 
analysis of radial intensity profiles, we have taken great 
care to remove uncertainties due to these imperfections. 

2.3. Image Reduction & Calibration 

We analyzed and reduced our data with the SCUBA 
User's Reduction Facility (SURF) (Jcnness & Lightfoot 
1997). The first-order effects of chopping had already been 
removed from the raw data, so our reduction entailed flat- 
fielding, extinction correction, removal of sky variations, 
and the coadding and rebinning of the maps. This pro- 
cess of data reduction is described in Paper I. The contour 
maps (with greyscale) for 850 and 450 /xm are in Figures 2- 
6. 

We calibrated our flux measurements through obser- 
vations of Uranus and AFGL 618, a primary and sec- 
ondary calibrator, respectively. We measured total volt- 
ages within 40" and 120" apertures. By dividing these val- 
ues into the expected flux from each of these sources^, we 
obtained the nightly calibration factors. For each night, 
we average the calibration factors and report this aver- 
age in Table 2. In February, we observed AFGL618 seven 
times and found significant variations (~40%) in our mea- 
surement of its flux at 450 /im. These variations are re- 
flected in the large standard deviation reported as uncer- 
tainties (i5(C|^°) and (5(C^|^)) in Table 2. To account for 
uncertainty in the flux calibration, we take the average 
of the calibration factors as determined from the obser- 
vations from all three nights in February and propagate 
the standard deviation as uncertainty in the flux measure- 
ment. This final average was used in the voltage-to-flux 
conversion for data obtained in February, 2000, and is re- 
ported in Table 2. 

2.4. Photometry 

We find the flux for each of our images by first mea- 
suring the total voltage within 40" and 120" apertures. 
Then, we use the calibration factors given in Table 2 to 
convert these values to flux measurements. The fluxes at 
450 and 850 /im arc given in Table 3. Statistical and cali- 
bration errors dominate our uncertainty (calculated as in 

^ We use the fluxes package in SURF to obtain the expected fluxe 
5460 = 173.37 Jy. The fluxes for AFGL 618 were SsbO = 4.56 Jy am 
page. These fluxes for AFGL 618 are comparable (within uncertainty) 



Paper I, Equation 2). In the modeling of these sources, we 
use the normalized intensity (/"°'''"(0)), so uncertainty in 
calibration has no effect on analysis of the radial profile. 

In our February 2000 observations, we had maps with 
a high signal-to- noise ratio, but the calibration errors ap- 
proach 100% at 450 /im (e.g., L1251B). In Table 3, we 
include the S/N ratio to give an indication of the quality 
of our maps. Additional photometry information is given 
in Table 4. The peak intensity centroids are given in Ta- 
ble 3. In addition, for IRAS 04016+2610 and 04361+2547, 
we give the fiux and centroid for neighboring condensa- 
tions. In Table 5, we give the calculated values for T(,o;, 

Lobsi and Lq})s/ Lgynm- 

2.5. Morphology 

We define an aspect ratio for each of the sources by di- 
viding the major and minor axis lengths as measured at 
the 2-cr contour. These values are reported in Table 5. 
For some cores (e.g., IRAS 03256+3055, 04016+2610 and 
04108+2803), the major and minor axes were not discern- 
able. In some cases, the cores are not ellipsoidal but, in- 
stead, have extensions in different directions. Also, the 
aspect ratio defined by the outer contours does not always 
reflect the morphology of the inner contours. The average 
aspect ratio was 1.3±0.4 with the most aspherical source 
being IRAS 04248+2612 with an aspect ratio of 2.3. 

For each modeled source, we determine the deconvolved 
source size (9dec) from the FWHM size of the beam {Omb) 
and source {Bsrc) intensity profiles: 6dec = V^src - ^'mb- 
The values for Odec/Smb are in Table 5; the reported un- 
certainties are based on the variability of 0mb iS0mb = 1- 
3"). In most cases (^^70%) , we resolve the emission, but 
the observed source size is sometimes comparable to the 
FWHM of the beam. For example, IRAS 04295+2251, 
with 9dec/0mb = 0.7, has an observed intensity profile 
that, at the FWHM, is only 1.3 times the FWHM of 
the beam (i.e., 9src = 20" and 0mb = 16"). In Fig- 
ure 1, we showed that the FWHM beamwidth does change 
slightly as the telescope cools. However, this mostly af- 
fects the 450 /im beam, and wc derive our deconvolved 
source sizes from the 850 /im profiles. Any source with 
6dec/Bmb > 1 is definitely resolved. Of course, those 
sources with 6dec/0mb = are clearly unresolved. For 
those sources with < Odec/dmb < 1, we have more closely 
analyzed their profiles. IRAS 04113+2758 was observed 
in November 1999 when the FWHM of the beam was 
quite variable (±3"), so we do not consider this a resolved 
source. IRAS 04016+2610, 04239+2436, and 04295+2251 
were all observed in February 2000, when the 850 fim beam 
was fairly stable. We consider each of these sources to be 
resolved. Further, we determine Odec/Smb from the FWHM 
of the intensity profiles, but our data is sensitive out to 3-5 
bcamwidths for most sources — i.e. we do resolve the outer 
envelope emission from these cores. As we will show in 
Section 3.1, the real discerning power of our observations 
is not in the inner regions but in the observations of emis- 
sion from the outer regions of the envelope. 

Also, several sources showed additional emission within 
the maps. For example, IRAS 04016+2610, 04108+2803, 

i for Uranus. In November 1999, these were: Ssbo = 65.63 Jy and 
5450 = 11-2 Jy as reported on the SCUBA secondary caUbrator web 
to the fluxes given, recently, by Jenness et al. (2002). 
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and 04361+2547 show sources of continuum emission that 
have comparable (or even greater) fluxes than the ob- 
served Class I protostar. Statistically, however, distinct, 
neighboring condensations of dust are not common in our 
sample — only about 20%. In Paper I, wc suggested Class 
I cores might be classified as such because of multiple 
sources within a beamwidth; in that sample, multiple 
sources were found in 9/16 of the observed Class and 
Class I cores. However, the occurrence of multiplicity is 
much less (3/16) in this sample of Class I cores. 

3. MODELS 

3.1. Power Law Models 

As in Paper III, we have modeled the dust emission 
from a subset of these Class I cores. We have selected 
the following sources for modeling based on their symmet- 
ric morphology and quality of data: CB230, L1251B, IRAS 
04166-F2706, IRAS 04169+2702, IRAS 04239+2436, IRAS 
04248+2612, IRAS 04295+2251, IRAS 04361+2547, and 
IRAS 04381+2540. We use a modification of the radia- 
tive transfer package by Egan, Leung, & Spagna (1988), 
CSdustS — a ray-tracing, 1-dimensional radiative transfer 
code. This program has several input parameters: lumi- 
nosity of the central source, opacity of the dust grains as a 
function of wavelength, strength of the interstellar radia- 
tion field (ISRF), and the density distribution of material 
surrounding the source. 

Primarily, we test power law density distributions, i.e., 

n(r) = n/ ; r € [r,, r^] (1) 

where n/ is the density at a fiducial radius, rt = 1000 
AU. The density has a lower limit of 1000 cm^"^, roughly 
simulating the ambient molecular cloud. We choose these 
models because they are simple and, a priori, make no as- 
sumptions about the theory incorporated by the collapse of 
material. In a later section, we do discuss the Shu collapse 
model (Shu 1977) because it has only one free parameter 
and is easily tested. Other theories of low-mass star forma- 
tion (e.g., Larson (1969), McLaughlin & Pudritz (1997); 
Foster & Chevalier (1993); Henriksen, Andre, & Bontemps 
(1997); Whitworth & Ward-Thompson (2001); etc.) often 
differ primarily in the value of p. 

For each of the model parameters, we have tested a 
reasonable range on two test cases: CB230 and IRAS 
04295+2251. Wc selected these sources as test cases be- 
cause of their high S/N maps and azimuthally symmetric 
morphology. The results from these tests are in Tables 6 
and 7. From these two test cases, we have also determined 
the uncertainty in p caused by the uncertainty in other 
parameters (sec Table 8). 

The luminosity of the central source (Lint) is determined 
by matching the modeled (Lmod) to the observed bolomet- 
ric luminosity {Lgbs)- At short wavelengths (A < 60/xrn), 
the cloud becomes optically thick, and the SED becomes 
more sensitive to opacity and geometry. Therefore, for 
Lmod and Lobs WC usc Only those flux points with A > 60 
/zm. For a given mass distribution (constrained by the 
intensity proflle) , we are able to determine the internal lu- 
minosity. Lint, required to produce Lots- This luminosity, 
Lint, may include both the luminosity of the protostar as 
well as accretion luminosity — basically, any source of lu- 
minosity within the inner radius. In most cases, Lm is 



less than Lots because the interstellar radiation flcld con- 
tributes a portion of the observed bolometric luminosity, 
~0.2 Lq. This is in accord with our results in Papers II 
and III. 

Wc usc dust opacities given by Osscnkopf & Henning 
(1994) for protostellar cores. They simulated the coagu- 
lation of dust grains (with an initial standard MRN size 
distribution) over 10^ yr and calculated the opacities for 
different initial gas densities. For the modeling of our sam- 
ple, we compare two "types" of dust from Table 1 of Os- 
senkopf & Henning. In Papers II & III, we have shown 
that the dust with thin ice mantles provides the best-fit 
models for prcprotostellar and Class sources. The opac- 
ities for this dust are given in the fifth column of Table 
1 in Ossenkopf & Henning (1994) and, hence, it is named 
"0H5" dust. In addition, we model sources with the opac- 
ities from column 2 of Table 1 (0H2 dust), representing 
dust grains without ice mantles. In the two test cases, we 
find that 0H5 dust fits the SED best, but the choice of 
dust type introduces a negligible source of uncertainty in 
p (see Table 8). 

We use the interstellar radiation field, given in Paper 11, 
based on Black (1994) and Draine (1978). For ultraviolet 
to far-infrared wavelengths, we introduce a factor, sjsrf, 
to scale the ISRF. For the two test cases, sisrf = 0.3 
fits the data best, in accord with the ISRF found to fit 
Class —1 and cores in Papers II & III. For this sam- 
ple, uncertainty in the ISRF contributes a large source of 
uncertainty in p: Sp ~ 0.4 (see Table 8). 

CSdustS requires the input of an inner and outer radius 
as boundary conditions to be used in solving the radia- 
tive transfer. We have developed some practical criteria 
in determining these radii. We set the inner radius to 100 
AU for all sources — this corresponds to 0'.'7 for the Taurus 
sources, 0'.'2 for CB230, and O'.'S for L1251B. The outer ra- 
dius should be, at least, the smn of the chop throw (120") 
and the radius of the model. For example, with CB230 
(at 450 pc), we used an outer radius of 100,000 AU which 
accounts for the 120" chop (54,000 AU) and for modeling 
out to 100" (45,000 AU). The outer radius for L1251B is 
72,000 AU; for Taurus sources, 33,000 AU. Wc do choose 
a lower limit for the density (n = 10^ cm~'^), but, for most 
cores, the outer radius lies within this limit. In Table 8, 
we show the effect of varying the inner and outer radii on 
the derived value of p. Changes in the inner radius cause 
a negligible change in p; uncertainty in the outer radius 
causes an uncertainty Sp = 0.3 in the modeling of CB230 
and 5p = 0.2 for IRAS 04295+2251. 

By assigning values to these various parameters and sim- 
ulating the transfer of radiation in the envelope, we calcu- 
late the temperature for the assumed density distribution. 
Then, as described in Paper III, we simulate the obser- 
vations including the effects of chopping and convolution 
of the data with the observed beam profile. The selec- 
tion of this beam profile was discussed in Section 2.2; we 
took great care in using the appropriate beam since it does 
vary with time. While a Gaussian beam will cause the de- 
rived value of p to be steeper by ^0.4, using an obseved 
beam profile from a different night causes negligible un- 
certainty (see Table 8). We attempted modeling of IRAS 
04302+2247, a source that was obviously unresolved. In 
this case, we were only able to determine that p > 3.5, 
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which, of course, wc expect in the modehng of a point 
source. 

We constrain the total amount of mass in the enve- 
lope by normalizing the modeled SED to the observed 850 
/im flux, and the luminosity of the internal source is con- 
strained by the observed luminosity, .Lobs- Then, we find 
the best-fit value for p, which we select by the lowest Y,y^ — 
the sum of the values for the 450 and 850 jim profiles 
(see Equation 4 in Paper III for the definition of Xa)- 
most cases, the best-fit model based on is identical 
to that selected by xlso (cf- L1251B, IRAS 04361, and 
IRAS 04381 in Table 9); because the uncertainty is much 
smaller for the 850 /im data, xiso much larger than xlso 
and, hence, dominates the sum. Even though the intensity 
profiles show Nyquist sampled data, we use, as in Papers 
II and III, only points spaced by a full beam to calculate 
X> • The derived value of p is virtually unaffected by the 
envelope mass or the luminosity of the internal source. We 
list the best-fit models in Table 9 along with those mod- 
els where p = ±0.1 and p = ±0.5. In Figures 7-15, we 
show the modeled and observed SEDs, radial profiles, and 
temperature distributions for each of our modeled sources. 

3.2. Top Hat Models 

We consider simple models of the density distribution 
in order to investigate the discerning capabilities of our 
methods — i.e. are we able to discern a power law distri- 
bution in the density from a top hat distribution? We 
define a top hat density distribution as an envelope model 
where the density is constant from the inner radius to some 
"break" radius; then, the density drops to 10"^ cm~^ in the 
region from this break radius to the outer radius. It should 
be most difficult to discern a power law from top hat dis- 
tribution with the source for which we derive the steepest 
power law, p = 2.3 for IRAS 04295+2251. In Figure 16, 
we show the modeled intensity profiles for three top hat 
distributions with breaks in the density at 5", 15", and 
25" (corresponding to 700, 2100, and 3500 AU). The max- 
imum density in the top hat functions is set so that the 
total mass is equal to that in the derived power law distri- 
bution. Also, in this figure, we show the best-fit power law 
model along with the observed intensity profile for IRAS 
04295+2251. This intensity profile shows a "bump" at 
about 5000 AU (or 35") that is coincident with additional 
emission within the field-of-view. We find the 5" and 15" 
top hat functions do not have enough emission at large 
radii to reproduce the 850 pxn profile. The 25" top hat 
has emission in the outer regions that is not seen in either 
the 850 or 450 /im observed profile. Likewise, the inten- 
sity profiles of other sources in our sample are not easily 
reproduced by simple top hat functions. Indeed, power 
law models, in contrast to the top hat functions, match 
the shape of the observed intensity profiles quite well. 

3.3. ShM Collapse Models 

In addition to modeling power law density distributions 
for our sample of nine cores, we have modeled these sources 
with the Shu inside-out collapse model (Shu 1977). This 
model begins with a singular isothermal sphere with a den- 
sity distribution that is proportional to r^^. Through 
some perturbation, collapse begins inside the cloud and 
proceeds outward at the sound speed (ae//)- As collapse 



ensues, the cloud's density distribution is described by an 
outer, static region with n oc and a region within the 
front of collapse (or infall radius, ri„/) whose density ap- 
proaches n oc r~^'^ (indicative of freefall). The beauty of 
this theory is that there is only one free parameter: time. 
As time progresses, the front of collapse (or infall radius, 
finf) propagates outward at a^ff- 

We determine a^f / from the observed linewidths in Ta- 
ble 10. Assuming a gas temperature of lOK, we are able 
to separate the turbulent and thermal components [aiurb 
and ath) of these linewidths and determine the effective 
sound speed [a^f j = V + ^turb) ■ ^'^^ e&ch source, ttef / 
is in Table 10. 

We list the parameters for our Shu collapse models in 
Table 11. For all modeled sources, the spectral energy dis- 
tribution is not fitted well by the Shu model. The Shu 
model overestimates all fluxes with A > 100 /zm, indicat- 
ing this model predicts too much mass within the cores. 
We may have overestimated the sound speed in these cores 
and, hence, also the mass (see Section 6.1). Also, varia- 
tions in dust opacity and geometry can greatly affect the 
SED (e.g., Moore & Doty 2002). Consequently, in choos- 
ing the best-fit Shu models, we give little regard to the 
values for Xsed- Instead, we focus on the nature of the 
intensity profiles wherein lies the strength of our analysis. 

The beam at 850 /im {9„ii, = 15") corresponds to 2100, 
4500, and 6800 AU for Taurus, L1251B, and CB230, re- 
spectively. For ah but two sources (IRAS 04169+2702 and 
04381+2540), the infaU radius of the best-fit Shu model 
falls within these radii. Basically, then, we are fitting 
straight power laws with an index oi p — 2. In Figure 17, 
we show the results of two models for IRAS 04295+2251. 
The best-fit model is with rinfaii = 500 AU, but we in- 
clude a model with rinfaii = 3000 AU. The larger infall 
radius causes the modeled profile to flatten. 

For IRAS 04169+2702, the best-fit infall radius is 
slightly larger than the JCMT beamwidth {rinfaii = 3000 
AU), and the best-fit Shu model has improved x^ values 
over the power law model {p = 1.5). In Figure 18, we 
show the results of this model. IRAS 04169+2702 was fit- 
ted with a power law of p = 1.5, so one might conclude 
that the Shu model should fit with an infall radius outside 
our modeling range. However, the Shu model predicts an 
inner region that approaches n oc r~^'^ . The density distri- 
bution just within the infall radius is actually significantly 
shallower than r~^'^, so the steep outer region is required 
to fit the data. 

Our basic conclusion for these Shu collapse models is 
that most sources in this sample exhibit no evidence for a 
break in the power law density distribution on the scales 
we can probe with SCUBA (~ 1000 - 8000 AU). 



3.4. Simple Analysis 

As in Paper I, we offer the results of a simple analysis 
involving the slope, m, of the outer part of the intensity 
profile (Adams 1991). Assuming the temperature distri- 
bution in the envelope to be described by a power law, 
Td(r) — Td{rf){r/rf)~'^ where r/ is a fiducial radius, one 
can determine the density power law index, which we call 
Pm, by the following relation: = m — q + 1 (see Equa- 
tions 5-8 in Paper I; assume q = 0.4). In Paper I, we 
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argued that this method is of dubious vaUdity for two rea- 
sons: failure of the Rayleigh- Jeans approximation and de- 
viation of the temperature profile from a power law due 
to heating from the ISRF and the central source. In or- 
der to compare with our more detailed modeling, we have 
performed such an analysis here to determine Pm- 

In Figures 7-15, we show the range over which we fit for 
m by the bold line on the x-axis of the 850 /xm intensity 
profile. The derived values oipm from the linear fits to the 
850 /im profiles arc in Table 12; we did not calculate m for 
the 450 /xm data. The average difference between pm and 
p is {pm —p)= 0.6 ± 0.3; i.e., the simple analysis tends to 
predict a much steeper density profile. This discrepancy is 
discussed in greater detail in Paper I (Section 4.2), but, in 
short, the failure of Rayleigh- Jeans approximation causes 
Pm to be steeper than the actual distribution; exclusion 
of the ISRF causes Pm to be shallower than the actual 
distribution. 

If we apply this correction to the sources analyzed in Pa- 
per I by this method, we find (pj^'') = 1.5 ±0.5 (compared 
to {p) = 1.6 ± 0.4 reported in this paper). 

Additionally, we apply the correction to the results of 
the large 1.3 mm survey by Motte & Andre (2001). At 
this wavelength, hv/k =11 K, so the Rayleigh- Jeans ap- 
proximation (where hu/kT -C 1) is better, but still not 
completely valid. Also, Motte & Andre chose values for 
q that ranged from —0.2 to 0.4 for their sample of Class 
and I cores. From Figures 7-15, it is clear that a neg- 
ative power law (indicative of external heating) does not 
describe the temperature profiles for these evolved proto- 
stars, and q = can be ruled out for all but one source, 
IRAS 04295-1-2251. Therefore, we assume q = 0.4 and cal- 
culate Pm from their values for m: {pm) — 1-8 ± 0.4, and, 
with our rough correction, {Pm"^) = 1.2dz0.5. However, be- 
cause the Rayleigh-Jeans approximation is less invalid at 
these longer wavelengths, this correction factor is probably 
too large and, hence, p^^ should be slightly larger. 

3.5. Disk Contribution 

We provide a simple analysis to determine the effect of 
the dust emission from a disk on our models. We assume 
the emission from a disk to be included only in the central 
beam of our observations. Mundy et al. (1996) found the 
disk for HL Tau to be <180 AU, and, since our beam sub- 
tends 2000-7000 AU, this assumption is reasonable. For 
the two test cases, CB230 and IRAS 04295+2210, we de- 
termine the effects of including a point source within the 
central beam that contributes one-half of the measured 
flux. In Figure 19, we show the results of adding this 
point source to the modeled intensity profile. The addi- 
tion of this point source causes the modeled intensity pro- 
file to steepen — hence, the best-fit intensity profile, when 
a point source is considered, can be matched by a shal- 
lower density profile. We find the best-fit density power 
law to be shallower by p = 0.5 when an unresolved com- 
ponent contributes one-half the fiux measured within the 
central JCMT beamwidth. Therefore, until the submil- 
limeter emission from disks is better constrained, we at- 
tribute an additional source of uncertainty in our calcula- 
tion of p: Sp = —0.5. However, this effect could be much 
more important for Class I cores, where the envelope is 
less substantial, than for Class cores. 



Observing submillimeter emission from disks requires 
submillimeter interferometers (e.g., SMA, ALMA), but we 
make some simple estimates of the disk emission for one 
source, IRAS 04361-1-2547. This source was observed in 
detail by Terebey et al. (1990, 1993) at near-infrared and 
millimeter wavelengths from which they suggest the exis- 
tence of a ~1000 AU disk or, at least, some dense compo- 
nent near the central position of this protostar. At A = 2.6 
mm, they report = 21 ± 7 mJy (1993) — observations 
at the Owens Valley Radio Observatory (OVRO) with a 
synthesized beam of 7'.' 7 x 6'.' 6. Further, they claim that 
these flux measurements are dominated by the emission 
from dust and not free-free emission. With our power law 
model for this source, we predict 5*,^ ^ 2 mJy (at A = 2.7 
mm); i.e., the observed flux does require the existence of 
an additional source of flux other than the envelope. In- 
deed, we consistently underestimate the long- wavelength 
flux (A > 1 mm) with our model of the envelope for all 
sources, further suggesting the existence of some \mre- 
solved component in these cores. 

We adopt a simple disk model (Butner et al. 1994) 
in order to estimate the 850 and 450 /um flux for IRAS 
04361+2547. The model consists of a central protostar 
surrounded by a disk whose surface density is described 
by a power law (E oc r~^'^). We assume a dust opacity of 
Ki, OC where /? = 1 and k;(850 /im)= 0.018 cm^gm"^. 
In Figure 20, we show the results of two calculations. The 
solid line represents a temperature profile in the disk that 
follows a power law, T oc r~° ''^. The dashed line shows 
a disk with T oc r~°'^; the outer edges of the disk are 
hotter and, thus, crudely simulate a flared disk. From the 
model of the "flared" disk, we set reasonable limits for the 
flux at 850 and 450 ^i^Jl■. S^{850 /im)< 0.4 Jy and S„(450 
/im)< 1.4 Jy. For each wavelength, these upper limits are 
60% of the observed flux in a 40" aperture. 

We have completed a similar analysis for IRAS 
04108+2803. Terebey et al. (1993) found 5^(2.6 mm)= 15 
mJy {9mb = 8'.'5 x 7'.'4) for this Class I source. Through 
the simple disk model, we set the following upper limits: 
5^(850 fim)< 0.24 Jy and S'^(450 fim)< 0.84 Jy. These 
upper limits correspond to 30% more than the observed 
850 /xm flux (within errors) and 74% of the observed 450 
/im flux. IRAS 04108+2803 is not resolved by our obser- 
vations, so it is conceivable that all of our observed flux 
is from a cental point source. Likewise, the observed 2.7 
mm flux for IRAS 04016+2G10 (6.8 mJy)(Hogerlieijde et 
al. 1997) sets an upper limit on the submillimeter fluxes 
originating from the presumed disk: 5^(850 /xm)< 0.15 Jy 
and 5*^(450 nm)< 0.73 Jy, 25% and 94% of the observed 
850 and 450 fluxes (40" aperture), respectively. The 
central component of IRAS 04016+2610 (L1489) was not 
resolved by our observations, so, again, it is possible that 
all of the submillimeter flux is from an unresolved cen- 
tral component. Indeed, since the nearby condensation to 
IRAS 04016+2610 lies partially within the 40" aperture, 
our reported flux values probably overestimate the flux of 
this Class I protostar. 

Chandler & Richer (2000) investigated the problem of 
disk contribution to the submillimeter continuum flux and 
concluded that the existence of a disk in Class sources, 
where the envelope is more substantial, has little effect on 
the derived value for p. Further, they claim the disk does 
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become a significant contributor in Class I objects. The 
envelope is less prominent in these cores, and a higher 
fraction of flux coming from the disk (i.e., by their for- 
malism, decreased fenv) causes a lower value for p. We 
conclude similarly, claiming that the inclusion of a disk 
causes the value for p to be decreased by 0.5; however, 
this effect could be significantly greater. Firmer conclu- 
sions await constraints from submillimeter interferometric 
observations. 

4. INDIVIDUAL SOURCES 
4.1. CB230 

CB230 (also IRAS 21169+6804) is the most luminous 
source wo observed. It is also the most distant at 450 pc 
(Launhardt & Henning 1997). Yun & Clemens (1994) ob- 
served a molecular outflow in the ^^CO 2-1 line for which 
the direction is shown on the 850 /xm contour map in Fig- 
ure 2. 

The submillimeter dust emission from CB230 is quite 
round with an aspect ratio of about 1. There does, how- 
ever, appear to be a slight extension in the direction per- 
pendicular to the outflow axis (oriented N-S). Due to the 
observed symmetry, we chose this source to model as a 
test case (see Table 6). The best-fit power law for CB230 
is p = 1.9. We find good fits for both the 850 and 450 
lj,m radial profiles, and the spectral energy distribution is 
matched well except for the 450 /im flux. This discrepancy 
is probably due to poor calibration at 450 /xm, but, since 
the focus of our modeling (i.e., the radial distribution) is 
independent of the absolute flux, we have not pursued this 
in great detail. 

4.2. IRAS 04295+2251 

IRAS 04295+2251 is the least luminous source we have 
modeled, with Lobs — 0.6 and Lint = 0.3 L©. It is 
clearly a Class I core by the bolometric temperature crite- 
rion with Tfcoi = 270 K, but it is Class by its submillime- 
ter luminosity Lobs/ Lgmm = 50. The core is resolved at 
850 /um, but its deconvolved size is the least of any of our 
modeled sources (^dec/^mh = 0.8). This source also has 
the highest derived value for p (2.3). We discuss the cor- 
relation between p and Odec/Omb in a later section. From 
the modeled density profile, we found the envelope mass 
to be 0.10 Mq. We chose this source as one of our test 
cases (see Table 7) because of its symmetric morphology 
(aspect ratio = 1) and high quality data (S/N=40 at 850 
jim). 

IRAS 04295+2251 is the only object in our sample that 
was detected by the NRAO VLA Sky Survey (NVSS) at 
21 cm (Condon et al. 1998). While it is not uncommon 
to find continuum radio emission (A > 3 cm) from Class 
I cores (Lucas et al. 2001), reports of strong emission at 
such a long wavelength are rare (cf. Skinner 1993). 

4.3. IRAS 04166+21/06 

IRAS 04166+2706 is a protostar associated with the 
B213 dark cloud in Taurus with a well-detected outfiow 
in the northeast to southwest direction (Bontemps et al. 
1996). The observations of IRAS04166+2706 were origi- 
nally presented in Paper I; but, since the core was classified 
as a Class I source (using the T^oi criterion), we model it 



here. The submillimeter emission is characterized by circu- 
larly symmetric contours except for an extension perpen- 
dicular (southeast) to the outfiow direction in the lowest 
contour interval (10% of the peak or 3(t; Paper I). We re- 
port the fluxes in a 120" aperture here: 5*850 = 1-9 ± 0.6 
Jy and 5450 = 12.9 ± 3.9 Jy. The bolometric luminosity is 
0.5 Lq and the bolometric temperature is 75 K. 

4.4. IRAS O4248+26I2 

IRAS 04248+2612 is about one arcminute north of 

the starless core Barnard 217 (B217) and is coincident 
with HH31 IRS2. Padgett et al. (1999) observed IRAS 
04248+2612 with HST/NICMOS in the near-infrared 
bands. Their images showed a clear bipolar structure, ori- 
ented NW-SE, that they claim is a reflection nebulosity, 
near-infrared light escaping by way of an outflow cavity. In 
addition, they observe a dark lane aligned perpendicular 
to this nebulosity. Our submillimeter maps (which exclude 
B217), show two sets of extensions: 1) on a small scale of 
~50", in the direction of the bipolar nebulosity (NW-SE) 
and 2) on a larger scale (^200"), aligned with the dark 
lanes observed by Padgett et al. (1999) in a NE-SW direc- 
tion. This large-scale structure is probably indicative of 
a filamentary-like distribution of mass in the cloud. The 
small-scale structure (within the central 30") is probably 
due to heating of the dust caused, in some way, by an 
outfiow. 

Myers et al. (1988) report no outflow activity in their 

^^CO 2-1 observations, for which the noise in one channel 
was ~0.2 K, but Moriarty-Schicven et al. (1992) find that 
the ^^CO 3-2 line observed at the central position shows 
low velocity line wings (Aw ~ 4.5 km s~^). Unfortunately, 
no molecular maps clarifying this outflow question have 
been published. 

Despite the asymmetric morphology of this source, we 
have attempted to model the dust emission. The best- 
fit power law is the lowest of any source, p = 0.8, and 
renders quite a good fit (see Figure 11). Of course, the 
extended emission is included in the azimuthal radial aver- 
age, so one might argue that this is the source of the mod- 
eled, flattened density distribution. However, we excluded 
the emission from the large-scale, NE-SW extensions and 
found that the best-fit power law distribution is still very 
low: p = 1.0. IRAS 04248+2612 is clearly different from 
all of our other sources, both in morphology and modeled 
density distribution. Its aspect ratio (see Table 5) is 2.3, 
the highest of any source. In Paper III, we showed that 
there does seem to be some correlation between aspect 
ratio and density distribution — i.e., those sources with a 
high aspect ratio tend to have flatter density distributions. 
It is possible that, in IRAS 04248+2612, we are observing 
the formation of a star under different physical processes 
and initial conditions than with the other cores. 

4.5. IRAS 03256+3055 

IRAS 03256+3055 is at 320 pc in the Perseus star- 
forming region. Unfortunately, this object has received 
little attention observationally, but we offer the first map 
of its dust emission (Figure 3). Ladd, Myers, & Goodman 
(1994) reported undetectable NH3 emission towards this 
object, but other reports of molecular line emission are 
found in the literature (Mardones et al. 1997, Gregersen 
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et al. 2000). We give its SED information in Table 4; 
this consists of a 60 ^ni flux and our submillimeter fluxes. 
Dent, Matthews, & Ward- Thompson (1998) report a very 
low 800 /zm flux — 0.125 Jy as compared to the 2.45 Jy 
we measure at 850 /xm. However, their observations were 
with a single bolometer with a beamsize of ~16" while the 
emission from IRAS 03256+3055 spans ~100". 

The luminosity of this object is 0.7 Lq, and, from its 
poorly sampled SED, we also derive a bolometric temper- 
ature of 16 K. This bolometric temperature is much lower 
than that reported by Jennings et al. (1987) (74 K) and, 
subsequently, Mardones et al. (1997), but they used IRAS 
upper limits on the flux at 12, 25, and 100 /im in their 
calculations. The envelope mass of IRAS 03256+3055 is 
1.7 M0 and is shared between, at least, three distinct dust 
condensations and significant extended emission. 

Nakamura & Li (2002) have conducted axisymmetric 
calculations concerning the formation of stellar groups via 
fragmentation. At least, in morphology, their results are 
strikingly similar to what we see in IRAS 03256+3055 — 
a ring- like structure with multiple condensations of dust. 
IRAS 03256+3055 provides an ideal laboratory in which 
to test the feasibility of such calculations as that presented 
by Nakamura & Li (2002) and, similarly, Boss (1993). 

4.6. IRAS 04016+2610, 04108+2803, 04361+2547, and 
04385+2250 

These four cores are located in Taurus and associated 
with L1489, L1495, TMRl, and TMCl. (respectively, 
for IRAS 04016+2610, 04108+2803, 04361+2547, and 
04385+2250). They are all Class I cores by their bolomet- 
ric temperatures (i.e., Tf^oi > 70 K), but the submillimeter 
luminosity of IRAS 04108+2803 [Lobs/Lsmm = HI) iden- 
tifies it as a Class core. 

Each of these cores has a neighboring condensation as 
shown in Figures 3 & 6. The fiux condensation associ- 
ated with IRAS 04016+2610 has been noted as a Class -1 
core by some (e.g., Hogerheijde & Sandell 2000), and we 
agree with this classification. Hogerheijde et al. (1998) 
observed an outflow associated with this Class I core that 
is not coincident with the flux condensation. 

IRAS 04361+2547 does exhibit a peak of submillimeter 
flux that is coincident with the outflow direction (outflow 
observed by Bontemps et al. 1996). Motte & Andre ob- 
served this source at 1.3 mm, but the condensation is not 
obvious at this longer wavelength indicative of the dust in 
this area being somewhat warmer than, for example, the 
Class —1 core associated with IRAS 04016+2610, which 
was detected by Motte & Andre (2000). Perhaps, the 
dust is heated directly by the outflow, or it is swept up 
and then heated by radiation from the central protostar 
that escapes by way of the outflow cavity. Regardless, 
IRAS 04361+2547 is probably not a Class -1 core. 

IRAS 04108+2803 and 04385+2250 arc both surrounded 
by significant extended emission, but, in contrast to 
04361+2547 and 04016+2610, this emission is not distinct 
from the central core. It could be that these cores have 
evolved by some asymmetrical collapse or that they are 
part of some larger star-forming complex. In the next 
section, we discuss IRAS 04385+2250 in the context of 
forming substellar objects. 



5. DISCUSSION 

5.1. Envelope Mass 

For the modeled sources, we have calculated the enve- 
lope mass by integrating the best-fit power law density 
distribution to an outer radius of 60" {M^nv per Equation 
6 in Paper III). Because the cloud is optically thin to sub- 
millimeter wavelength radiation, the following expression 
is is often used to estimate the envelope mass (Hildebrand 
1983): 

ISO J'^L' 

where S„ is the flux density at 850 /xm in a 120" beam 
(Table 3), i?^ is the Planck function, Kj, is the opacity 
per gram of gas and dust at 850 /xm {ki, = 1.8 x 10~^ 
cm^^gm^^ for 0H5 dust), Tiso is the isothermal dust tem- 
perature, and D is the distance to the source. Of course, 
the dust is not isothermal, but we use our models to esti- 
mate the isothermal temperature (Tiso) that makes Equa- 
tion 2 agree with the masses derived from the more de- 
tailed models (see Equation 7, Paper III). We use the re- 
sulting average Tiso (16±4 K) for the modeled sources to 
calculate the isothermal envelope mass (Af,!*" ) for the un- 
modelcd sources. The values for M,,nv and Af^fj" arc in 
Table 10. In Figure 21, we show a histogram of these cal- 
culated masses for the Class I sources in this paper and the 
Class cores in Papers I and IV. The envelopes of Class 
I cores are clearly less massive than Class sources. Only 
CB230, of the Class I cores, has M > 0.5 M©, and it is 
very close to the Class boundary with T^oi = 74 K. The 
remaining cores have M < 0.5 Mq and are in Taurus. 

In Table 10, we also give the FWHM linewidth, associ- 
ated reference, and the calculated virial mass. The virial 
mass, is corrected for the best-fit power law density 

distribution as in Paper III (Equations 8 & 9). For those 
sources not modeled, we adopt the average p (1.6) to cal- 
culate the virial mass. To be consistent with our dust mass 
calculations, we adopt an aperture of Gap = 120", which 
corresponds, in Taurus, to 16,800 AU. We plot the ratio of 
virial to dust mass {M^r / ^dust) against Mdust and T{,o; 
in Figure 22. For Tboj > 70 K, these two calculations of 
the envelope mass begin to diverge. Reasonable variations 
in the dust opacity can introduce a factor of ten uncer- 
tainty (Ossenkopf & Henning 1994), but these variations 
cannot account for values of M^r /Mdust ^ 60. One sce- 
nario is that the molecular tracers used, mostly N2H+, do 
not trace the material surrounding Class I cores. Caselli 
et al. (2002) mapped two Class I cores coincident with our 
sample: IRAS 04016+2610 and IRAS 04248+2612 (listed 
as L1498 and B217, respectively, by Caselli et al.). In both 
cases, the nearby preprotostellar core is observed to have 
strong emission in these tracers, but the tracer shows no 
peak on the Class I core. Simply, observations of N2H+ 
do not trace the mass in these Class I cores. 

5.2. Luminosity and Stellar Mass 

As discussed in the previous section, the masses for the 
envelopes of these cores are very low. Likewise, the ob- 
served luminosities are very low. For those cores in Tau- 
rus, we find a range of bolometric luminosities from 0.2 
to 4.6 Lq with an average of 1.4+1.3 L©. These lumi- 
nosities are upper limits on the luminosity of the central 
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protostellar system because of the ISRF contribution. For 
the most luminous sources, L1251B and CB230, the ISRF 
contributes negligibly to the bolometric luminosity, but, 
for those cores in Taurus, the ISRF accounts for a consid- 
erable fraction of Lots ■ Wc estimate the total contribution 
of the ISRF to the observed luminosity to be about 0.2 
L0. 

Based on the standard picture of a Class I protostar 
in which the observed luminosity results from accretion, 
wc can calculate the mass of this protostar by a simple 
relationship (e.g.. Stabler 1994): 

GM 

Here, M* is the stellar mass. Lace is the luminosity due to 
accretion, G is the gravitational constant, R^, is the stellar 
radius, and M is the accretion rate given as M = noCi^/G, 
where qt is the isothermal sound speed and Uo is a nu- 
merical constant of order unity (Stabler 1994). Assuming 
a minimum isothermal sound speed of 0.2 km s~^, we cal- 
culate M = 2 X 10"^ Mgyr-^ Then assuming i?* = 3 R©, 
an overestimate of the stellar radius based on the models 
of Stabler (1988), we can calculate an upper limit on the 
stellar mass from this relationship. For the least luminous 
of the cores in this sample (IRAS 04385+2250; Lots = 0.2 
Lq), the mass, based on Equation 3, is A'/* = 0.01 Mq 
or, since this in the substellar realm, = 10 Mj„p. The 
envelope mass for this core is 0.04 Mq, and, assuming a 
star formation efSciency of 25% (Wilking et al. 1989), the 
final mass of this object would be Af* = 0.02 M© (or 20 
Mjup) in about 10** yr (assuming constant M = 2 x 10~^ 
Mgyr-i). 

This analysis holds many shortcomings and uncertain- 
ties. The mass accretion rate (A^) is highly uncertain, 
possibly episodic, and perhaps 1 — 2 orders of magnitude 
smaller than that used in this calculation (e.g., MuzeroUe 
et al. 1998). Such an accretion rate would increase Af* by 
a factor of 10 — 100. On the other hand, the observed lumi- 
nosity is probably some combination of the disk accretion 
luminosity and luminosity of the central protostar. Indeed, 
MuzeroUe et al. (1998) find, in some cases, the accretion 
luminosities contribute ~ 10% to the total observed lumi- 
nosity. If only 10% of the total luminosity results from 
accretion (as estimated by MuzeroUe et al. 1998), the cal- 
culated stellar mass is a factor of 10 smaller. 

We can also estimate the stellar mass by other means. If 
we assume an age (r), we can place IRAS 04385+2250 on 
calculated isochrones. The models of Siess et al. (2000), 
with T = 10^ yr, predict ~ 0.1 Mq. Burrows et al. 
(2001) give ~ 0.2 Mq for r = 10*^ yr. These models 
probably overestimate M^, however, because they assume 
no accretion and have not considered the initial condi- 
tions governed by infall and collapse of the protostellar 
envelope. Also, the age of embedded Class I cores has 
been estimated to be ~ 10^ yr (e.g., Myers et al. 1987, 
Kenyon et al. 1990); most stellar evolution models begin 
at r = 10^ yr and are highly uncertain for ages less than 
lO'^ yr. Siess et al. (1999) included accretion {M = 10"'^ 
MQyr"""^) in their calculations of pre-main sequence evo- 
lutionary tracks, and, for t — 4.5 x 10^ yr, their models 
would indicate Af* < 0.1 Mq; unfortunately, these models 
do not include objects with masses lower than 0.1 Mq as 
required to accomodate IRAS 04385+2250. Additionally, 



the application of these models to Class I cores is quite 
precarious; these models are intended for T Tauri stars. 

Observations of young, forming (Class O/I) brown 
dwarfs have yet to be firmly established, but the exis- 
tence of older substellar objects with disks (e.g. Natta & 
Testi 2001) suggest that brown dwarfs do form by similar 
mechanisms as T Tauri stars. If so, we, indeed, probably 
have observed brown dwarfs within their natal envelopes. 
Then, it should not be surprising that, in this sample of 
low-luminosity Class I cores, a fraction arc substellar. In- 
deed, Kenyon & Hartmann (1995) concluded similarly, i.e. 
the "luminosity problem" for Class I objects is solved if 
such objects are substellar. Of course, these authors also 
present another viable solution, that the accretion radius 
is larger in these young objects than we expect. 

5.3. Classification: Lots, Tboii and Lamm 

Currently, the class system (Class — 1 — III) implies a 
sequence of stages that are occupied at some point by 
all evolving protostars (e.g., Andre & Montmerle 1994). 
Hence, statistical arguments have been constructed by 
comparing the numbers of sources in each stage to estab- 
lish the relative time spent in each epoch (e.g., Andre & 
Montmerle 1994). As discussed previously, three criteria 
for classification have been put forth: unir (Lada et al. 
1987), Tboi (Myers & Ladd 1993), and Lobs/Lsmm (Andre 
et al. 1993). These values are given in Table 5 for our 
sample of sources. 

For the selection of the SED points from which these 
values were calculated (see Table 4), we required that 
the aperture be large enough to include all emission: 
9mb > 20". This restriction causes the exclusion of the 
near-infrared (NIR) fluxes, but we also calculate T^oi with 
the NIR data because the bolometric temperature is most 
affected by the NIR data. The NIR flux, no matter how 
weak, causes this mean frequency to be skewed to higher 
values. Indeed, if wc do not include the NIR flux in our 
calculation of T(,o;, 35% of our Class I sources become Class 

by the T^oi criterion. Of course. Class I sources are iden- 
tified by their NIR emission and are even called, by some, 
NIR protostars. However, detection of the NIR light is 
often dependent on orientation and, hence, could feasibly 
be obscured in an otherwise Class I object. 

An additional complication in distinguishing between 
Class I and cores is found in the luminosity scheme de- 
veloped by Andre et al. (1993). They define a Class I 
source as one for which Lobs/Lsmm > 200, assmning that, 
as a protostar evolves, the envelope should be cleared away 
and, hence, this ratio should increase. Originally, this in- 
dicator was based on measurements within a small aper- 
ture (20") because of nearby sources. However, for our 
isolated cores, it is clear that such a small aperture does 
not sample the total submillimeter flux emitted by the ob- 
ject. Therefore, we use only large aperture measurements 
(> 20") in the calculation of these values (Table 5); for 
our submillimeter fluxes, we used the 40" aperture fluxes 
in the calculation of Lobs/Lsmm so that nearby condensa- 
tions would not contribute to the submillimeter flux. 

In Figure 23, we show a plot of Tboi versus Lobs/Lsmm 
that includes all sources from Paper I and this paper. Class 

1 cores are characterized by Tboi > '^OK (Myers & Ladd 
1993, Chen et al. 1995) or Lobs/Lsmm > 200 (Andre et 



10 



Young et al. 



al. 1993). We have subdivided the plot into four regions 
based on these boundaries. Class I cores should be found 
only in the upper right-hand portion of the plot, and Class 

cores should be in the lower left-hand portion. This is 
clearly not the case. While there is a noticeable correla- 
tion between the two evolutionary indicators, the divisions 
between the classes are not so obvious. Indeed, all of the 
cores in this paper's sample are Class I by the Ttoi cri- 
teria (except for L1251B and IRAS 03256+3055); half of 
these Class I cores, however, arc Class by the Lobs/ Lgmm 
criterion. The problem arises in that the T^oi observ- 
able is greatly affected by geometry whereas Lobs/Lsmm 
is not so affected. Visser, Richer, & Chandler (2002) also 
found such a discrepancy between these classification crite- 
ria. They find that the submillimeter luminosity criterion 
{Lobs/ Lsmm > 200) renders far fewer Class I cores than 
does the bolometric temperature. Further, they infer life- 
times for the evolutionary stages based on the numbers of 
Class and I cores. In order to create such statistical ar- 
guments concerning the lifetimes of Class — l-III objects, 
one must recognize that different numbers of Class or I 
objects will be found depending on the system of classifi- 
cation employed. 

5.4. Density Structure 

The primary purpose of our modeling is to derive the ac- 
tual distribution of matter in the envelopes of these form- 
ing protostars. For the Class I cores in this paper, we find 
the density power law exponent {p) = 1.6 ± 0.4 (mean ± 
the standard deviation of the sample); however, the pres- 
ence of a disk or point-like component could cause p to 
decrease by as much as 0.5. The median for this sam- 
ple of Class I cores is p = 1.8. In Paper III, we reported 
(p) = 1.6 ± 0.3 for the Class sources. For the Class 
and I objects in Papers III and IV, we find (p) = 1.6 ±0.4, 
and the median is p = 1.8. If we leave out those Class or 

1 sources with aspect ratios > 1.5, a total of four sources 
from Papers III and IV, we find (p) = 1.8 ± 0.2, and the 
median is also p = 1.8. We include these modeled sources 
in a histogram (Figure 24) both separately and together. 
The value for p is clearly peaked between 1.5-2.0, seem- 
ingly inconsistent with models that predict very shallow 
distributions — e.g., McLaughlin & Pudritz (1997). How- 
ever, these models should not be excluded altogether, for 
we do sec evidence, in some sources, of a shallower profile, 
and the inclusion of a disk allows for significantly shallower 
density profiles. 

In Figure 25, for all modeled sources (from Paper III 
and IV), we show the aspect ratio as a function of p. 
A clear correlation exists between these two values such 
that the more aspherical sources exhibit shallower pro- 
files. Of course, the azimuthal averaging of an aspheri- 
cal object could result in a shallow intensity profile (and, 
hence, shallower density profile), but we claim this is a 
small effect. In Section 4.4, we discussed the modeling of 
IRAS 04248-1-2612 and how the derived density power law 
changes negligibly with exclusion of the extended emission 
in the NE-SW extensions. We also found this to be true 
in aspherical, Class sources — see Paper III. 

As in Paper III, we find little evidence for variation of 
dust properties within the envelope. We compare the nor- 
maUzed intensity ratios (liio/IL/ISio!!^) Figure 26 and 



conclude that, for most sources, the observed variations in 
this ratio are well-fitted by including the effects of external 
heating by the ISRF and convolution of the actual beam. 
For L1251B, IRAS 04361+2547, and IRAS 04381+2540, 
we are unable to simultaneously match the observed 850 
and 450 /um intensity profiles. The dust opacity may vary 
with radius in these sources, but it could also be that we 
have not appropriately modeled the temperature in these 
cores, perhaps, due to asymmetric heating, externally, by 
the ISRF or, internally, by a disk+protostar or binary cen- 
tral source. 

While power law density distributions do not have a 

physical size, we define an observational size from the de- 
convolved FWHM of our modeled sources (Table 5). In 
Figure 27, we show the derived value for p plotted against 
the effective size {9dec/ Smb) for the modeled Class and I 
cores (from Papers III & IV). Additionally, the solid line 
shows a scries of models for 0.5 < p < 2.5 for which the 
deconvolved model size is measured with the 15" beam. 
The models show that steeper power laws arc clearly less 
well-resolved in agreement with the observed trend. 

The trend in Figure 27, while not surprising, is quite 
striking; smaller sources exhibit a larger value for p. We 
stress this point for three reasons. First, Terebey, Chan- 
dler, & Andre (1993) addressed this resolution issue — 
that steeper power law density distributions are less well- 
resolved. However, one must be very careful when ana- 
lyzing the intensity profiles of marginally resolved sources 
because uncertainty in the beamshape can be ~ 1 — 3" 
(see Section 3.2). Indeed, there are examples of marginally 
resolved sources for which the authors have derived very 
steep density power laws {p > 2.5), density distributions 
that have little theoretical basis. 

Also, we propose an intriguing possibility implying this 
trend is not entirely systematic and might be indicative 
of some physical quality. Perhaps, we are observing only 
the outer envelope of these less well-resolved cores. That 
is, these "smaller'' sources are not small because they are 
more evolved and have collapsed over time, but they are 
small because they began with an initial preprotostellar 
core that was lower in mass and, hence, has a smaller outer 
radius than the other cores in our sample. In the Shu col- 
lapse model, we would expect a steeper density profile in 
the outer, static envelope, and it might be that, indeed, 
we are seeing just that. 

A third scenario is also possible. Perhaps, the "smaller" 
cores, which we interpret to have steeper density profiles, 
have a greater disk contribution. This point source of emis- 
sion, as shown in a previous section, results in the deriva- 
tion of a steeper power law distribution in the density. In 
this case, the envelope mass distribution could be consid- 
erably shallower that we predict here. Indeed, if the disk 
contribution is more significant as the core size decreases, 
the trend for the envelope density distribution shown in 
Figure 27 could be shifted such all cores harbor an enve- 
lope described by 1.0 < p < 1.5, a possibility that would 
suggest the validity of the theory presented by McLaughlin 
& Pudritz (1997). 

Wc compare these results to those recently presented by 
J0rgensen, Schoier, and van Dishoeck (2002) (hereafter, 
JSD) and find that we reach similar conclusions. These 
authors model two sources coincident with our sample: 



Tracing the Mass during Low-Mass Star Formation, IV 



11 



IRAS 04361+2547 and IRAS 04381+2540, which they caU 
TMRl and TMCl, respectively. Our derived values of the 
density power law exponent for these two sources are con- 
sistent with those of JSD. However, their modeling differs 
from ours in the assigning of two parameters: the ISRF 
and outer radius; additionally, they do not simulate the 
effects of chopping. Exclusion of the ISRF should cause 
JSD to derive a shallower power law exponent to fit the 
observed data. Also, JSD adopt much smaller outer radii 
than we do in this paper- in some cases, a factor of four 
smaller. By assigning a small outer radius and excluding 
the effects of chopping, these authors should derive a some- 
what steeper power law for the density distribution than 
that presented here. The combination of these two effects 
(i.e., exclusion of the ISRF and truncated outer radius) 
seem to cancel in their analysis, rendering similar results 
to our own. 

5.5. Caveats 

As presented in the previous sections, we have discussed 
several caveats that hinder direct interpretation of the 
dust continuum emission: beam convolution of our model, 
complex temperature distributions due to external heat- 
ing, different assumed dust opacities, contribution to the 
observed flux from a disk component, and the inner and 
outer radii assumed in the model. We do not analyze the 
effects of outflows because, for most of our sources, the out- 
flows are either very weak or undetected (e.g., Bontemps 
et al. 1997). Our modeling is also limited by its one- 
dimensional nature, but we point out objects for which 
three-dimensional modeling will prove fruitful (e.g., IRAS 
04248+2612). Our assertion that asymmetric sources ex- 
hibit shallower density distributions suggests that these 
objects are forming in truly different environments. How- 
ever, this issue will not be fully resolved until models of 
radiative transfer for three dimensions are completed. 

6. CONCLUSIONS 
Our conclusions are as follows: 

1) We find that 3/16 (20%) of the sources in 
this sample have distinct, neighboring condensations 
(IRAS 03256+3055, 04016+2610 and 04361+2547). 
These condensations could be the result of com- 
plex structure (03256+3055), neighboring preprotostel- 
lar cores (04016+2610), or the heating due to outflows 
(04361+2547). 

2) We have fit power law models for the density dis- 
tribution in the outer envelopes of nine Class I cores and 
find (p) = 1.6 ± 0.4 and a median of p = 1.8. The exis- 
tence of a disk or point- like component could cause p to be 
much shallower by about 0.5. Table 12 has a summary 
of parameters for the best-fit models. 

3) We test the Shu inside-out collapse model and find 
that, in most cases, this model provides a good fit for the 
intensity profile but with an infall radius that is within 
the central resolution element of our observations. Most 
sources in this sample exhibit no evidence for a break in 
the power law density distribution on the scales we can 
probe with SCUBA (- 1000 - 8000 AU). 

4) We compare our analysis with a simpler analysis 
put forth by Adams (1991) and as used in Paper I. This 
method assumes a temperature power law (with index 



q = 0.4) and hv/kT <^ 1. These assumptions are invalid 
for these low luminosity sources, and this simple analysis 
overestimates p by ~0.6. 

5) In modeling the density distribution, we have deter- 
mined several sources of uncertainty: possible existence of 
a disk {6{p) = -0.5), variations in the ISRF {5{p) = ±0.4), 
and uncertainty in the outer radius {5{p) = ±0.3). Con- 
tributions by the disk to the observed submillimeter flux 
will be addressed by such instruments as SMA and ALMA 
while the ISRF may be constrained by future observations 
in the mid- to far- infrared (e.g., SIRTF). We find no ev- 
idence for a truncated outer radius in models completed 
over a wide range of radii. 

6) For our modeled sources, we have calculated the 
mass of the envelope surrounding the protostar: {M^nv) = 
6.6 ± 2.4 M© for the two most luminous sources and 
{Menv) = 0.2 ± 0.1 Mq for those modeled sources in Tau- 
rus. We find that these masses are clearly less than those 
found for Class cores (Paper HI) and, further, contradict 
the virial masses predicted by N2H+ linewidths. Molecu- 
lar line tracers such as N2H+ simply do not trace the mass 
in these Class I cores. 

7) The standard classification criteria, Tboi and 
Lobs/Lsmm, are sometimes inconsistent for Class and 
I cores. The mode of classification is important for stud- 
ies that attempt to estimate core lifetimes based on the 
Class and Class I source counts (e.g. Visser, Richer, & 
Chandler 2002). 

8) Sources with higher aspect ratios are better fitted 
with shallower density power laws. While this effect is 
partially due to our one-dimensional models, we propose 
that, indeed, these sources are forming in physically dif- 
ferent initial conditions. Three-dimensional modeling will 
help to unravel this question. 

9) We note a trend between deconvolved FWHM source 
size and p — smaller cores have higher values for p. This 
trend is expected since steeper density power law distri- 
butions produce less well-resolved intensity profiles. It is 
also possible that we are observing the outer portions of 
the envelope in these less well-resolved sources. 

10) We have also discussed several soiirces individually. 
In particular, IRAS 03256+3055, with its unique morphol- 
ogy, is an ideal object on which to test theories of frag- 
mentation in the formation of low-mass protostars. Also, 
we note the possibility, through some simple calculations, 
that IRAS 04385+2550 is a young, forming substellar ob- 
ject with a mass of ~ 10 Mj„p. 
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Fig. 1. — The observed JCMT beam profiles were obtained through observations of Uranus or AFGL618 during November 1999 and 
February 2000, respectively. For each month, we show the beams from one representative night as measured at different times (given in 
Hawaii- Aleutian Standard Time). The solid lines are Gaussian functions with a FWHM of 14" and 8" for 850 and 450 /im, respectively. 
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Fig. 2. Maps of emission at 850 and 450 fim. For all the maps, the lowest contour is at tlie 2-a level. The contours increase by intervals of 
ten percent of the peak for CB230 (10-cr and 2-ct for 850 ^m and 450 /im) and L1251B (6-(t and 2-a for 850 fini and 450 ^m). The contours of 
IRAS 04302+2247 are shown in increments of twenty percent (lO-cr) for the 850 /im map and ten percent {4-(t) for the 450 /im map. Central 
positions are in Table 1. Outflow directions are labeled with bold lines; references for outflows are as follows; CB230, Yun & Clemens, 1994.; 
L1251B, Sato et al., 1994; IRAS 04302+2247, Bontemps, et al. 1996. 
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FlH. 3. Maps of emission at 850 and 450 fim. For all the maps, the lowest contour is at the 2-cr level. The 850 ;tm eontours inerease by 
intervals often pereent of the peak for IRAS 03256+3055 (l-cr), IRAS 04016+2610 (4-cr), and IRAS 04108+2803 (2-0-). Contours for IRAS 
04016+2610 at 450 ^im are drawn at levels of 10 pereent (2-ct). The 450 fim maps for IRAS 03256+3055 and IRAS 04108+2803 have eontours 
at 1-a intervals. Central positions are in Table 1. Outflow directions are labeled with bold lines; references for outflows are as follows; IRAS 
04016+2610, Hogerheijde, et al. 1998. 
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Fig. 4. — Maps of emission at 850 and 450 fim. For all the maps, the lowest contour is at the 2-a level. The 850 fini contours increase by ten 
percent of the peak (IRAS 04113+2758, 6-cr; IRAS 04169+2702, 3-a; IRAS 04239+2436, 2-a). The 450 contours for IRAS 04113+2758 
and IRAS 04169+2702 increase by twenty percent (4- and 6-cr, respectively) and by ten percent {2-a) for IRAS 04239+2436. Central positions 
are in Table 1. Outflow directions are labeled with bold lines; references for outflows are as follows: IRAS 04169+2702, Bontemps et al., 
1996.; IRAS 04239+2436, Saito et al., 2001. IRAS 04113+2758 shows no detectable outflow activity (Saito et al. 2001). 
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Fig. 5. — Maps of emission at 850 and 450 /im. For all the maps, the lowest contour is at the 2-a level. For IRAS 04248+2612 and IRAS 
04264+2433, the 850 and 450 fim contours increase by l-cr intervals. IRAS 04295+2251 has contours that increase by ten percent of the peak 
(2-<t) at 850 /im and twenty percent (3-(7) at 450 /im. Central positions are in Table 1. 
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Fig. 6. Maps of emission at 850 and 450 fini. For all the maps, the lowest contour is at the 2-cr level. The 850 and 450 fim contours 
increase by ten percent of the peak for IRAS 04361+2547 and IRAS 04381+2540 corresponding to 5- and 2-ct for IRAS 04361+2547 (850 and 
450 ^m, respectively) and 1-ct for IRAS 04381+2540 (both maps). The contours for IRAS 04385+2550 arc drawn at the 1-cr level. Central 
positions are in Table 1. Outflow directions are labeled with bold lines; references for outflows are as follows: IRAS 04361+2547 and IRAS 
04381+2540, Bontemps et al. 1996. 
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Fig. 7. — Model of CB230. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy distribution 
is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and observed radial 
profiles — the best-fit model is a solid dark line. Also shown in these plots arc the beam profiles used in the beam convolution (solid light 
line). The range of fit for m (in the calculation of pm) is shown by the bold line on the x-axis. 
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Fig. 8. — Model of IRAS 04166+2706. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy 
distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution 
(solid light line). The range of fit for m (in the calculation of pm) is shown by the bold line on the x-axis. 
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Fig. 9. — Model of IRAS 04169+2702. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy 
distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution. 
The range of fit for m (in the calculation of pm ) is shown by the bold line on the x-axis. 
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Fig. 10. — Model of IRAS 04239+2436. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy 
distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution 
(solid light line). The range of fit for m (in the calculation of pm) is shown by the bold line on the x-axis. 
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Fig. 11. — Model of IRAS 04248+2612. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy 
distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution 
(solid light line). The range of fit for m (in the calculation of pm) is shown by the bold line on the x-axis. 
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Fig. 12. — Model of IRAS 04295+2251. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy 
distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution. 
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Fig. 13. — Model of IRAS 04361+2547. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy 
distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution 
(solid light line). The range of fit for m (in the calculation of pm) is shown by the bold line on the x-axis. 
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Fig. 14. — Model of IRAS 04381+2540. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy 
distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution 
(solid light line). The range of fit for m (in the calculation of pm) is shown by the bold line on the x-axis. 
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Fig. 15. — Model of L1251B. In the upper left panel, we show the modeled temperature profile. The modeled spectral energy distribution 
is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and observed radial 
profiles — the best-fit model is a solid dark line. Also shown in these plots arc the beam profiles used in the beam convolution (solid light 
line). The range of fit for m (in the calculation of pm) is shown by the bold line on the x-axis. 
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Fig. 16. — Wc show the modeled intensity profiles for the best-fit power law of IRAS 04295+2251 (p = 2.3, solid line) along with a series 
of top hat functions. The top hat functions (dashed and dotted lines) have steps at 5", 15", and 25" where the density drops to 10^ cm~^. 
The bottom panels also have the observed beam profile shown by the light, solid line. 
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Fig. 17.— Shu collapse models for IRAS 04295+2251 with two different infall radii, 500 (solid) and 3000 AU (dashed). In the upper left 
panel, we show the modeled density profile. The modeled spcetral energy distribution is marked by crosses in the upper right panel, the 
observed fluxes are circles with errorbars. The two bottom plots have the modeled and observed radial profiles — the best-fit model is a solid 
dark line. Also shown in these plots are the beam profiles used in the beam convolution (solid light line). 
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Fig. 18. — Shu collapse model for IRAS 04169+2702. In the upper left panel, we show the modeled density profile. The modeled spectral 
energy distribution is marked by crosses in the upper right panel, the observed fluxes are circles. The two bottom plots have the modeled and 
observed radial profiles — the best-fit model is a solid dark line. Also shown in these plots are the beam profiles used in the beam convolution 
(solid light line). 
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Fig. 19. — The observed radial intensity profiles are shown as error bars for CB230 and IRAS 04285+2251. The best-fit models, assuming 
all emission originates in only the envelope, are p = 1.9 and p = 2.3 for CB230 and IRAS 04295-1-2251, respectively. The solid line represents a 
model of the envelope that has a shallower density profile by 0.5 (p = 1.4 and p = 1.8). The dashed line results from adding a point source to 
this modeled envelope emission. This point source can cause the modeled intensity profile from the envelope+disk to be considerably steeper, 
and, hence, the derived density power law for the envelope is shallower. We assume moderate fluxes for the disk and find the effect to be a 
decrease in p by 0.5, but the eflfect could possibly be much more significant. 



32 



Young et al. 



1 - 







^-1 



o 



-Z 



-3 - 



-4 



"I 1 1 1 r 

IRAS 04361+2247 




O.lJy < S^(850/xm) < 0.4Jy 



0.2Jy < S„(450Mm) < 1.4Jy 



\ \ I 







Log A (/xm) 



Fig. 20. — The modeled SED of a disk surrounding a protostar. The surface density profile of the dust within this disk is a power law 
(S oc r~^-'^). The solid line shows a model where the temperature power law is steep (T oc r"''-''^^) while the dashed line model has a shallower 
profile (T oc r^"'^). The latter is intended to simulate the effects of flaring. The crosses are observed 2.6 mm flux density from Terebey et al. 
(1993) and the submillimeter fluxes in this paper (for 40" aperture). The ranges for Si, at 850 and 450 /im as derived from the two models 
are given — these ranges are reasonable upper and lower limits on the dust emission from the disk at these wavelengths. 
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Fig. 21. — The histogram of envelope masses for Class I and Class cores from this paper and Paper III. We show the mass derived from 
the modeled density profile for the modeled sources and the isothermal mass (assuming 100:1 gas-to-dust ratio) as discussed in the text. All 
Taurus cores in this sample have Menu < 0.5 Mq. Those sources with M^nv > 0.5 Mq have T^oi < 75 K. 
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Fig. 22. — We plot the ratio of virial to dust mass {Myir/^dust) against M^ust and T},oi. Class I cores are represented by squares; Class 
cores (from Paper III) are circles. 
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Fig. 23. — The bolometric temperature (Tj,o() is plotted against LobsjJ-'sm.m- The dashed lines represent the boundary between Class 
and I. Class I objects have T^,„l > 70 K (Myers & Ladd 1983) and L„,,^/Lsmm > 200 (Andre et al. 1993). 
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Fig. 24. — All modeled sources from Paper III and this paper axe represented in the top histogram. Class and I (from Paper III and this 
paper) sources are shown in the bottom and middle histograms, respectively. 
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Fig. 25. — The aspect ratio (given in Table 5) is on the horizontal axis and the modeled density power law exponent (p) is on the vertical 
axis. More asphcrical objects are best modeled with shallow density profiles — an effect that is possibly real and not just an artifact of our 
one-dimensional model. 
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Fig. 26. — The error bars represent the observed ratio of 450 to 850 /im flux. The solid line represents the modeled value. 
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Fig. 27. — The derived value for p is plotted against the ratio of deconvolved source size and beam size (at FWHM). The solid line represents 
models of the dust emission with 0.5 < p < 2.5. We note an obvious trend of increasing p for the less well-resolved sources. This trend is 
possibly an effect due to the cores compact nature, but we also note the possibility that we are observing the outer parts of the envelope in 
these "smaller" cores. 
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Table 1 
Observed Sources 



Source a (B1950.0) 5 (B1950.0) a (J2000.0) 5 (J2000.0) Date Source 
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(° 
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"\ 
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(h 
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S 


I 


(° ' 


) 




cn vi r o nmcnt 


CB230 


21 


16 


55, 


.0 


68 


04 


52 


21 


17 


40, 


,0 


68 


17 


31.9 


11/19/99 


IRAS 21169+6804 


L1251B 


22 


37 


40, 


.8 


74 


55 


50 


22 


38 


47, 


,2 


75 


11 


28.8 


11/19/99 


IRAS 22376+7455 


IRAS 03256+3055 


03 


25 


39, 


.2 


30 


55 


20 


03 


28 


44, 


,5 


31 


05 


39.7 


02/02/00 


Per 


IRAS 04016+2610 


04 


01 


40, 


.0 


26 


10 


49 


04 


04 


42, 


,5 


26 


18 


57.9 


02/03/00 


L1489 


IRAS 04108+2803 


04 


10 


49, 


.3 


28 


03 


57 


04 


13 


54, 


,9 


28 


11 


30.5 


02/03/00 


L1495 


IRAS 04113+2758 


04 


11 


20 


.9 


27 


58 


30 


04 


14 


26, 


,4 


28 


06 


01.4 


11/19/99 


L1495 


IRAS 04166+2706'' 


04 
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37 


.8 
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29 


04 


19 
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,5 


27 


13 


39.7 


08/30/98 




IRAS 04169+2702 


04 
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.0 
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02 


52 


04 
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,6 


27 


10 


01.7 


02/03/00 


B213 


IRAS 04239+2436 


04 
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.5 


24 


36 


54 


04 


26 
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.4 


24 


43 


35.9 


02/02/00 


L1524 


IRAS 04248+2612 


04 


24 


53, 


.2 


26 


12 


39 


04 


27 


57, 


,2 


26 


19 


16.9 


02/01/00 


HH31 IRS/B217 


IRAS 04264+2433 


04 


26 


28, 


.1 


24 


33 


24 


04 


29 


30, 


,0 


24 


39 


55.6 


02/03/00 


L1524 


IRAS 04295+2251 


04 


29 


32, 


.2 


22 


51 


11 


04 


32 
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,1 


22 


57 


30.3 


02/01/00 


L1536 


IRAS 04302+2247 


04 


30 


16, 


.6 


22 


47 


05 


04 


33 


16, 


,4 


22 


53 


21.3 


02/02/00 


L1536 


IRAS 04361+2547 


04 


36 
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.4 
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47 


27 


04 


39 
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53 
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TMRl 


IRAS 04381+2540 


04 
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.6 


25 


40 


48 
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41 


11, 


,6 


25 


46 


32.1 
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IRAS 04385+2550 


04 


38 


34, 


.2 


25 


50 


43 


04 


41 


38, 


.4 


25 


56 


25.3 


02/03/00 


TMCl 



^CB230 is at a distance of 450 pc (Launhardt & Henning 1997), and L1251B is at 300 pc (Kun & Prusti 1993). 
Perseus sources (a = 03 h) are located at a distance of 320 pc (de Zeeuw et al. 1999) while Taurus sources (a = 04 

h) are at 140 pc (Elias 1978). 

''The observation of IRAS 04166+2706 was reported in Paper I, but we reanalyze it here. 



Table 2 



Observing and Calibration Summary 



Date 


T850 


7450 


'^CSO 


r'SSO 
'^40 


'-'120 


'-'40 


r'450 
'-'120 










Jy/V- 


Jy/V^ 


Jy/V- 


Jy/v*^ 


November 
















11/19/99 


0.27 (0.08) 


1.24(0.09) 


0.061(0.006) 


0.94 (0.08) 


0.78 (0.03) 


2.73 (0.78) 


2.02 (0.46) 


February 
















02/01/00 


0.23 (0.03) 


1.23(0.18) 


0.048(0.007) 


0.83 (0.04) 


0.68 (0.03) 


1.87(0.63) 


1.44(0.95) 


02/02/00 


0.17 (0.03) 


0.79(0.27) 


0.044(0.002) 


0.84 (0.01) 


0.67 (0.01) 


2.54(0.26) 


1.71(0.31) 


02/03/00 


0.21 (0.02) 


1.02(0.16) 


0.050(0.006) 


0.83 (0.03) 


0.69 (0.03) 


2.88(0.09) 


2.36(0.06) 










February Average 














0.83(0.03) 


0.68(0.03) 


2.35(0.61) 


1.78(0.71) 



^Calibration Factor for a 40" diameter aperture 
Calibration Factor for a 120" diameter aperture 
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Table 3 



Observed Flux Densities of Sources 



Source 


Centroid 


S,. (Jy) 


(Jy) 


(Jy) 


S. (Jy) 


"450/850'' 


"450/850*' 


S/N'^ 




(Aa",A5") 


850 iJ.m''- 


850 nm^ 


450 iJ.m'^ 


450 nm^ 








CB230 


(-8,2) 


2.74(0.33) 


3.03(0.28) 


5.13(1.76) 


7.76(2.11) 


1.0 


1.5 


110 


L1251B 


(4,3) 


3.66(1.50) 


6.69(2.70) 


15.04(14.25) 


24.31(7.07) 


2.2 


2.0 


70 


IRAS 03256+3055 




0.45(0.05) 


2.45(0.29) 


1.90(1.93) 


8.08(3.89) 


2.3 


1.9 


30 


IRAS 04016+2610 


(10,0) 


0.59(0.03) 


1.92(0.10) 


4.23(1.38) 


11.13(6.65) 


3.1 


2.8 


40 




(74,16) 


0.78(0.03) 




3.53(1.38) 




2.4 






IRAS 04108+2803 


(-3,-9) 


0.17(0.02) 


0.61(0.06) 


1.13(0.92) 


5.38(2.31) 


2.9 


3.4 


14 


IRAS 04113+2758 


(0,-1) 


1.06(0.12) 


1.65(0.15) 


1.86(0.63) 


4.36(1.27) 


0.9 


1.5 


180 


IRAS 04166+2706 


(1,-4) 


1.08(0.06) 


1.9(1.0) 


4.2(1.0) 


12.9(6.0) 


2.1 


3.0 


18 


IRAS 04169+2702 


(4,2) 


1.14(0.05) 


2.53(0.12) 


6.09(1.88) 


11.75(5.05) 


2.6 


2.4 


60 


IRAS 04239+2436 


(2,-2) 


0.42(0.05) 


0.78(0.10) 


2.14(2.26) 


3.72(1.79) 


2.6 


2.5 


50 


IRAS 04248+2612 


(-2,1) 


0.56(0.10) 


1.95(0.33) 


2.96(3.03) 


8.97(3.93) 


2.6 


2.4 


20 


IRAS 04264+2433 


(3,-3) 


0.13(0.01) 


0.37(0.04) 


0.63(0.47) 


1.32(0.57) 


2.5 


2.0 


10 


IRAS 04295+2251 


(1,-6) 


0.42(0.02) 


0.78(0.05) 


2.66(1.01) 


5.34(2.34) 


2.9 


3.0 


40 


IRAS 04302+2247 


(5,1) 


0.57(0.08) 


0.62(0.08) 


2.09(2.50) 


1.88(0.91) 


2.0 


1.8 


90 


IRAS 04361+2547 


(5,4) 


0.64(0.08) 


1.21(0.11) 


2.35(0.80) 


4.55(1.32) 


2.1 


2.1 


60 




(63,-57) 


0.48(0.08) 




1.69(0.80) 




2.0 






IRAS 01.381+25.10 


(15.6) 


0.56(0.01) 


1.35(0.11) 


2.82(1.13) 


6.11(2.80) 


2.0 


2.5 


10 


IRAS ()l:-!8.")+2."').")n 


((iA.T) 


0.08(0.01) 


().2<)(().01) 


0.11(0.11) 


0.37(0. 10) 


0.8 


0. 1 


10 



''In a 40" aperture 

''In a 120" aperture 

"S/N is given for the 850 /im maps. 
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Table 4 

Spectral Energy Distributions of Sources 



Source 


A (//m) 


Sjy (Jy) 


^mb (") 


Ref. 
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300x45 


10^ 




25 


3.88(0.39) 


300x45 


10^ 




60 


7.49(0.75) 


90x300 


10^ 




100 


10.86(0.11) 


180x300 


10^ 




450 


5.38(2.31) 


120 


11^ 




800 


0.085(0.033) 


16.8 


2 




850 


0.61(0.06) 


120 


11^ 




1100 


<0.1 


18.5 


2 




1300 


<0.02 


60 


14 


IRAS 04113+2758 


2.2 


0.32 (0.05) 


10.7 


3^ 




2.2 


0.66(0.02) 


9 


13 




12 


2.03 (0.28) 


300x45 


10^ 
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Table 4 — Continued 



Source 


A (/im) 


Su (Jy) 


^mh (") 


Ref. 




60 


12.75 (0.13) 


90x300 


10" 




100 


15.82 (2.21) 


180x300 


10" 




450 


4.36(1.27) 


120 


11" 




800 


0.98(0.12) 


16.8 


2 




850 


1.65(0.15) 


120 


11" 




1100 


0.461(0.053) 


18.5 


2 




1300 


0.410(0.040) 


11 


14 




1300 


0.750( ■ ■ ■ ) 


60 


14" 




3400 


0.052(0.006) 


8 


17 


IRAS 04166+2706 


1.6 


0.00010(0.00002) 


10 


4 




2.2 


0.00019(0.00009) 


10 


4 




12 


0.07(0.007) 


300x45 


10'* 




25 


0.58(0.058) 


300x45 


10* 




60 


5.9(0.59) 


90x300 


10* 




100 


9.5(0.95) 


180x300 


10* 




450 


12.9(6.0) 


120 


11* 




850 


1.9(1.0) 


120 


11* 


IRAS 04169+2702 


2.2 


0.02(0.003) 


10.7 


3* 




2.2 


0.02(0.004) 


9 


13 




12 


0.75(0.05) 


300x45 


10* 




25 


5.21(0.31) 


300x45 


10* 




60 


17.00(1.7) 


90x300 


10* 




100 


17.46(2.4) 


180x300 


10* 




450 


11.75(5.05) 


120 


11" 




800 


0.75(0.18) 


16.8 


2 




800 


0.52(0.069) 


13.5 


13 




800 


0.73(0.051) 


16.8 


13 




850 


2.53(0.12) 


120 


11" 




1100 


0.280(0.025) 


18.5 


13 




1100 


0.281(0.053) 


18.5 


2 




1300 


0.190(0.009) 


11 


14 




1300 


0.730( ■ ■ ■ ) 


60 


14" 




3400 


<0.026 


8 


17 


IRAS 04239+2436 


2.2 


0.07 (0.01) 


10.7 


3" 




12 


1.71(0.1) 


300x45 


10" 




25 


6.98(0.42) 


300x45 


10" 




60 


15.24(0.15) 


90x300 


10" 




100 


15.88(2.22) 


180x300 


10" 




450 


3.72(1.79) 


120 


11" 




800 


0.333(0.043) 


16.8 


2 




800 


0.23(0.025) 


13.5 


13 




850 


0.78(0.10) 


120 


11" 




1100 


0.077(0.023) 


18.5 


13 




1100 


0.114(0.021) 


18.5 


2 




1300 


0.080(0.01) 


11 


14 




1300 


0.170 


60 


14* 




3400 


<0.024 


8 


17 


IRAS 04248+2612 


2.2 


0.04(0.006) 


10.7 


3* 




2.2 


0.04() 


3.8 


16 




25 


1.33(0.13) 


300x45 


10* 




60 


4.62(0.46) 


90x300 


10* 




100 


9.26(0.93) 


180x300 


10* 




450 


8.97(3.93) 


120 


11* 




800 


0.252(0.046) 


16.8 


2 




800 


0.12(0.033) 


13.5 


13 




800 


0.21(0.039) 


16.8 


13 




850 


1.95(0.33) 


120 


11" 




1100 


0.100(0.019) 


18.5 


13 




1100 


0.099(0.015) 


18.5 


2 




1100 


0.15(0.045) 




18 




1300 


0.060(0.007) 


11 


14 




1300 


0.450 


60 


14" 


IRAS 04264+2433 


2.2 


0.026(0.001) 


9 


13" 




12 


0.39(0.02) 


300x45 


10" 




25 


2.94(0.18) 


300x45 


10" 




60 


5.21(0.52) 


90x300 


10" 




450 


1.32(0.57) 


120 


11" 




850 


0.37(0.04) 


120 


11* 
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Table 4 — Continued 



Source 







^mh \ } 


Ref. 


1300 


0.031(0.008) 


11 


14 


1300 


0.031( ■ ■ ■ ) 


60 


14a 


2.2 


0.12 (0.02) 


10.7 


3^ 


12 


0.60(0.11) 


300x45 


10^ 


25 


1.82(0.18) 


300x45 


10^ 


60 


3.53(0.35) 


90x300 


lO''' 


100 


7.41(1.04) 


180x300 


10=^ 


450 


5.34(2.34) 


120 


11* 


800 


0.241(0.050) 


16.8 


2 


850 


0.78(0.05) 


120 


IP 


1100 


0.094(0.018) 


18.5 


2 


1300 


0.115(0.010) 


11 


14 


1300 


0.115( •■■ ) 


60 


14a 


2.2 


0.03(0.005) 


10.7 


3* 


2.2 


0.025(0.001) 


9 


13 


2.2 


0.025(...) 


3.8 


16 


25 


0.44(0.06) 


300x45 


10* 


60 


6.40(0.13) 


90x300 


10* 


100 


9.43(0.19) 


180x300 


10* 


450 


1.88(0.91) 


120 


11* 


800 


0.342(0.057) 


16.8 


2 


850 


0.62(0.08) 


120 


11=^ 


1100 


0.149(0.019) 


18.5 


2 


1300 


0.180(0.010) 


11 


14 


1300 


0.180( •■■ ) 


60 


14a 


2.2 


0.04(0.006) 


10.7 


3^ 


2.2 


0.04(0.001) 


9 


13 


12 


1.81(0.11) 


300x45 


10^ 


25 


18.87(1.13) 


300x45 


10^ 


60 


44.75(6.27) 


90x300 


10=^ 


100 


35.43(4.96) 


180x300 


10* 


450 


4.55(1.32) 


120 


11* 


800 


0.634(0.067) 


16.8 


2 


850 


1.21(0.11) 


120 


11* 


1100 


0.188(0.027) 


18.5 


2 


1300 


0.110(0.008) 


11 


14 


1300 


0.440( ■ ■ ■ ) 


60 


14* 


2.2 


0.01 (0.002) 


10.7 


3* 


12 


0.44(0.06) 


300x45 


10* 


25 


2.69(0.27) 


300x45 


10* 


60 


10.29(1.0) 


90x300 


10* 


100 


13.91(1.4) 


180x300 


10* 


450 


6.41(2.80) 


120 


11* 


450 


12.5() 




15 


800 


0.289(0.053) 


16.8 


2 


850 


1.35(0.11) 


120 


11* 


850 


2.60() 




15 


1100 


0.116(0.013) 


18.5 


2 


1300 


0.070(0.009) 


11 


14 


1300 


0.300( ■ ■ ■ ) 


60 


14* 


12 


0.54(0.03) 


300x45 


10* 


25 


1.55(0.16) 


300x45 


10* 


60 


2.88(0.29) 


90x300 


10* 


450 


0.37(0.16) 


120 


11* 


850 


0.29(0.04) 


120 


11* 


1300 


0.030(0.005) 


11 


14 


1300 


0.030( ••• ) 


60 


14* 



IRAS 04295+2251 



IRAS 04302+2247 



IRAS 04361+2547 



IRAS 04381+2540 



IRAS 04385+2550 



Flux value used in the calculation of Tjoj and Lob^ 



References. — 1. Myers et al. 1987;2. Moriarty-Scheiven et al. 

1994; 3. Tamura et al. 1991; 4. Kcnyon ct al. 1990; 5. Launhardt 
et al. 1997; 6. Huard et al. 1999; 7.Morcira ct al. 1997; 8. Yun 
et al. 1995; 9.Kun & Prusti 1993; 10. IRAS PSC, 1988; 11. this 
paper; 12. Launhardt & Hcnning, 1997; 13. Barsony & Kenyon, 
1992; 14. Motte & Andre, 2000; 15. Hogerheijde & Sandell, 2000; 
16. Padgett et al., 1999; 17. Saito et al., 2001; 18. Dent et al., 1998. 
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Table 5 
Observed Source Properties 



Source 




Tbol 


Lobs 


Lsmm/ Lobs 


Aspect 


6dec/6mb 




K 


K 


Le 




a 
b 






CB230 


47 


74 


11.5 


100 


1.2+0.2 


1.3+0.3 


L1251B'' 


64 




10.8 


100 


1.7+0.2 


1.7+0.4 


IRAS 03256+3055^ 


16 




0.7 






6.9+0.9 


IRAS 04016+2610 


92 


167 


4.6 


1000 




0.7+0.3 


IRAS 04108+2803 


98 


179 


1.0 


500 




0.0+0.3 


IRAS 04113+2758 


82 


274 


2.0 


333 




0.7+0.3 


IRAS 04166+2706 


50 


75 


0.5 


50 


1.2+0.2 


1.4+0.3 


IRAS 04169+2702 


83 


133 


1.4 


111 


1.2+0.2 


1.9+0.2 


IRAS 04239+2436 


103 


194 


1.8 


500 


1.4+0.3 


0.8+0.3 


IRAS 04248+2612 


49 


136 


0.9 


142 


2.3+0.4 


3.0+0.3 


IRAS 04264+2433 


109 


186 


0.5 


500 




0.0+0.3 


IRAS 04295+2251 


91 


270 


0.6 


125 


1.1+0.2 


0.8+0.2 


IRAS 04302+2247 


53 


122 


0.5 


111 




0.0+0.3 


IRAS 04361+2547 


98 


141 


3.6 


1000 


1.0+0.2 


1.1+0.2 


IRAS 04381+2540 


78 


126 


0.9 


142 


1.2+0.2 


1.7+0.2 


IRAS 04385+2550'' 


122 




0.2 


500 




0.0+0.3 



*Near-infrared data is not included in the calculation of this value. 
''No near-infrared emission has been detected towards these objects. 
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Table 6 
Test Case, CB230 



Source p Uf n To Unt SISRF xlbO X850 XSED 

(xlO^cm-^) (AU) (xlO^AU) Lq 



ISRF 



1.9 


22 


100 


1 


0H5 


6.7 


0.3 


0.4 


28 


91 


2.0 


22.5 


100 


1 


0H5 


5.8 


0.6 


0.9 


31 


105 


2.2 


25 


100 


1 


0H5 


4.5 


1 


0.6 


53 


113 


2.2 


< 25 


100 


1 


0H5 




3 








Inner Radius 




















2.2 


30 


50 


1 


0H5 


4.5 


1 


0.6 


53 


113 


2.2 


25 


100 


1 


0H5 


4.5 


1 


0.6 


53 


113 


2.2 


25 


200 


1 


0H5 


4.75 


1 


0.5 


49 


120 


2.2 


25 


400 


1 


0H5 


6.75 


1 


0.3 


43 


117 


2.2 


25 


800 


1 


0H5 


6.8 


1 


1 


38 


200 


Outer Radius 




















1.9 


35 


100 


0.25 


0H5 


4.8 


1 


1.6 


33 


32 


2.1 


33 


100 


0.5 


0H5 


4.7 


1 


1.9 


50 


28 


2.1 


27 


100 


0.75 


0H5 


4.4 


1 


2.9 


48 


88 


2.2 


25 


100 


1 


0H5 


4.5 


1 


0.6 


53 


113 


2.2 


0.5 


100 


2 


0H5 


6.0 


1 


33 


560 


591 


Dust Type'' 




















2.2 


25 


100 


1 


OH5 


4.5 


1 


0.6 


53 


113 


2.2 




100 




OH2 




1 








Beam'' 




















2.2 


25 


100 


1 


0H5 


4.5 


1 


2 


20 


113 



''We were unable to fit the luminosity or fiducial density of CB230 with 0H2 dust, but the 

intensity profiles were fit well with p — 2.2. 

''We use the beam profile observed from the February observations (CB230 was observed in 
November). 
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Table 7 
Test Case, IRAS 04295+2251 

Source p nj r\ ro Lint sjsrf xho xiso Xsed 

(xlO^cm-^) (AU) (xlO^ylf/) Lq 



ISRF 



2.3 


2 


100 


3.3 


0H5 


0.25 


0.3 


0.3 


2 


12 


2.5 


1.5 


100 


3.3 


0H5 


0.19 


0.6 


0.3 


4 


12 


2.7 


0.8 


100 


3.3 


0H5 


0.1 


1 


0.3 


5 


18 


2.7 


< 0.8 


100 


3.3 


0H5 


0.1 


3 








Inner Radius 




















2.7 


1 


50 


3.3 


0H5 


0.06 


1 


0.3 


5 


20 


2.7 


0.8 


100 


3.3 


0H5 


0.1 


1 


0.4 


4 


23 


2.6 


1 


200 


3.3 


0H5 


0.1 


1 


0.6 


3 


21 


2.6 


1 


400 


3.3 


0H5 


0.4 


1 


0.6 


4 


25 


Outer Radius 




















2.5 


3.6 


100 


1.1 


0H5 


0.19 


1 


0.4 


2 


16 


2.5 


2.2 


100 


2.2 


0H5 


0.12 


1 


0.9 


2 


16 


2.3 


0.8 


100 


3.3 


0H5 


0.1 


1 


3.3 


3.5 


21 


2.3 




100 


6.6 


0H5 




1 








Dust Type 




















2.3 


0.9 


100 


3.3 


OH2 


0.05 


1 


6 


5 


24 


2.3 


0.8 


100 


3.3 


OH5 


0.1 


1 


3.3 


3.5 


21 


Beam'' 




















2.7 


0.8 


100 


3.3 


0H5 


0.1 


1 


2 


5 


19 



''We use the beam profile observed from the November observations (IRAS 04295+2251 was 
observed in February). 



Table 8 

Modeling Uncertainty Analysis 



Variable 


Range 


Ap 


ALint 


ISRF 


0.3-3'* 


±0.3 


2.2 




0.3-3'^ 


±0.4 


0.15 


Inner Radius 


50-400AU'' 


< 0.1 


2.3 




50-400AU'' 


< 0.1 


0.3 


Outer Radius 


25,000-200,000AU'' 


±0.3 


1.2 




11,000-66,000AU'' 


±0.2 


0.1 


Dust Type 


OH2-OH5 


< 0.1 




Disk 


20-80% of total flux 


-0.5 




Beam 


Nov99-Feb00 


< 0.1 






Uncertainty:'' 
Uncertainty:'^ 


+U.4 
-0.7 
+0.4 
-0.7 


3.4L0 
0.4Lq 



"For CB230. Outer radii range corresponds to 55 — 
440". The uncertainty is simply the quadratic sum of the 
contributing factors and, since each is not entirely inde- 
pendent, this uncertainty is probably an overestimate. 

''For IRAS 04295+2251. Outer radii range corre- 
sponds to 80 - 430". 
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Table 9 
Power Law Models'' 



SourcG 


n 






A450 


A860 


2 

X.SED 


C o rnniGiit s 


CB230 


2.0 


30 


6.6 


1 


45 


37 






1.9 


22 


6.7 


0.4 


28 


91 


Best-fit 




1.8 


18 


7.2 


1 


37 


119 






2.4 


55 


6.8 


10 


242 


169 






1.4 


6.0 


12 


34 


455 


251 




IRAS04166+2706 


1.7 


3.2 


0.4 


0.2 


3 


64 






1.6 


3.0 


0.4 


0.2 


2 


72 


Best-fit 




1.5 


2.4 


0.4 


0.3 


3 


92 






2.1 


6.6 


0.2 


2 


16 


57 






1.1 


1.0 


0.4 


6 


26 


28 




IRAS 04169+2702 


1.6 


3.5 


0.73 


11 


50 


35 






1.5 


2.8 


0.77 


10 


57 


39 


Best-fit 




1.4 


2.5 


0.80 


10 


71 


40 






2.0 


7.0 


0.55 


50 


94 


18 






1.0 


0.8 


0.9 


155 


308 


57 




IRAS 04239+2436 


2.0 


1.0 


1.0 


1.7 


1.6 


94 






1.9 


1.0 


1.0 


0.8 


0.7 


18 






1.8 


0.9 


0.9 


0.1 


2 


6 


Best-fit 




1.7 


0.75 


0.9 


0.2 


5 


12 






2.3 


1.5 


0.75 


5 


9 


3 






1.3 


0.3 


1.7 


18 


66 


150 




IRAS 04248+2612 


0.9 


0.5 


0.5 


0.7 


1 


11 






0.8 


0.35 


0.6 


0.5 


2 


12 


Best-fit 




0.7 


0.25 


0.95 


1.4 


6 


13 






1.3 


1.5 


0.4 


7 


20 


5 






0.3 


0.09 


1.0 


15 


67 


22 






l.O'^ 


1.0 


0.5 


1 


3 


8 




IRAS 04295+2251 


2.4 


2.0 


0.25 


0.3 


4 


11 






2.3 


2.0 


0.25 


0.3 


2 


12 


Best-fit 




2.2 


2.0 


0.25 


0.9 


2 


12 






2.8 


2.0 


0.2 


2 


13 


7 






1.8 


0.9 


0.4 


9 


14 


29 




IRAS 04361+2547 


1.9 


1.9 


2.2 


9 


17 


16 






1.8 


1.8 


2.2 


6 


14 


17 


Best-fit 




1.7 


1.5 


2.2 


4 


16 


20 






2.3 


3.0 


1.5 


17 


49 


6 






1.3 


0.5 


3.0 


8 


125 


33 




IRAS 04381+2540 


1.7 


2.0 


0.53 


2 


17 


4 






1.6 


1.8 


0.53 


4 


9 


4 


Best-fit 




1.5 


1.4 


0.53 


7 


4 


6 






1.4 


1.0 


0.6 


12 


1 


8 






1.3 


0.8 


0.7 


22 


1 


9 






2.0 


3.2 


0.6 


1 


54 


1 






1.0 


0.35 


1.1 


78 


42 


15 




L1251B 


1.4 


8.3 


13 


22 


183 


646 






1.5 


12 


10 


41 


111 


73 


Best-fit 




1.6 


16 


10 


61 


172 


82 






2.0 


44 


7.5 


138 


1050 


1261 






1.0 


2.3 


14.0 


36 


2033 


570 





^Wc set an inner radius of lOOAU for all sources. The outer radius is assigned as follows: 
CB230, 100,000AU; L1251B, 72,000AU; Taurus sources, 33,000AU. The ISRF is 0.3 times 
the standard as discussed in the text. 

''We excluded the extended emission for this model. See text for a more detailed discus- 
sion. 
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Table 10 



Derived Masses & Velocities 



Source 


AvFWHM 

km s~^ 


km s ^ 


Mq 




T 

ISO 

K 


Av 
Rcf. 


CB230 


0.8 


0.39 


15 


4.60 


11 


2 


L1251B 


1.00 


0.46 


18 


4.95 


14 


1 


IRAS 03256+3055 


0.34 


0.23 


6.1 


1.74 


16 


1 


IRAS 04016+2610 


0.23 


0.21 


2.4 


0.26 


16 


3 


IRAS 04108+2803 


0.22 


0.20 


2.3 


0.08 


16 


4 


IRAS 04113+2758 


0.22= 


0.20 


2.3 


0.22 


16 


4 


IRAS 04166+2706 


0.34 


0.23 


2.3 


0.34 


13 


1 


IRAS 04169+2702 


0.29 


0.22 


2.4 


0.37 


15 


1 


IRAS 04239+2436 


0.26^= 


0.21 


2.4 


0.08 


19 


4 


IRAS 04248+2612 


0.34^= 


0.23 


4.3 


0.14 


23 


4 


IRAS 04264+2433 


0.26^^ 


0.21 


2.4 


0.05 


16 


4 


IRAS 04295+2251 


0.3'^ 


0.22 


1.1 


0.10 


16 


4 


IRAS 04302+2247 


0.3^= 


0.22 


2.5 


0.08 


16 


4 


IRAS 04361+2547 


0.36^= 


0.24 


2.5 


0.16 


16 


4 


IRAS 04381+2540 


0.36^= 


0.24 


2.8 


0.18 


15 


4 


IRAS 04385+2550 


0.27= 


0.22 


2.4 


0.04 


16 


4 



^Mv"'' is calculated by Equation 8 in Paper III. 

''The envelope mass, Me„t,, is based on the best-fit density distribution (per 
Equation 6 in Paper III) from which an isothermal temperature is derived. The av- 
erage isothermal temperature (16K) is then used for the other sources to calculate 
the envelope mass, M*^" , via the method derived by Hildebrand (1983) (Equation 
4, Paper I). 

=These velocity linewidths are for the region in which these cores reside. How- 
ever, their positions are often several beamwidths from the IRAS source (~3-10'). 

References. — 1. Mardones et al. 1997 (N2H+); 2. Wang et al. 1995 (C^^O); 
3. Fuller & Myers 1993 (HC3N); 4. CaseUi et al. 2002 (N2H+). 
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Table 11 
Shu Collapse Models 





'^infall 
ATT 


km. s 




v2 

A450 


v2 

A8B0 


X-SED 


1 .ilTTl m OTl'i'G 
WlllllldlljO 


CB230 


1000 


0.39 


6.7 


2 


21 


4600 






2000 


0.39 


6.7 


0.8 


18 


2089 


Best-fit 




3000 


0.39 


6.7 


0.7 


78 


1251 






4000 


0.39 


6.7 


1 


170 


880 




IRAS 04166+2706 


1000 


0.23 


0.4 


0.3 


3 


43 






2000 


0.23 


0.4 


0.2 


0.8 


36 


Best-fit 




3000 


0.23 


0.4 


0.8 


4 


33 




IRAS 04169+2702 


1000 


0.22 


0.8 


23 


58 


92 






2000 


0.22 


0.8 


6 


29 


72 






3000 


0.22 


0.8 


6 


29 


64 


Best-fit 




4000 


0.22 


0.8 


12 


49 


60 




IRAS 04239+2436 


500 


0.21 


0.9 


0.5 


1 


682 


Best-fit 




1000 


0.21 


0.9 


0.3 


5 


488 






2000 


0.21 


0.9 


3 


20 


328 




IRAS 04248+2612 


10,000 


0.23 


0.6 


6 


13 


146 


Best-fit 




15,000 


0.23 


0.6 


7 


17 


140 




IRAS 04295+2251 


500 


0.22 


0.25 


3 


12 


928 


Best-fit 




1000 


0.22 


0.25 


10 


23 


759 




IRAS 04361+2547 


500 


0.24 


2.2 


8 


14 


359 






1000 


0.24 


2.2 


4 


12 


655 


Best-fit 




1500 


0.24 


2.2 


2 


18 


546 




IRAS 04381+2540 


3000 


0.24 


0.5 


10 


3 


187 


Best-fit 




4000 


0.24 


0.5 


14 


2 


152 






5000 


0.24 


0.5 


17 


4 


131 




L1251B 


2000 


0.46 


10 


52 


161 


2316 






3000 


0.46 


10 


36 


23 


1989 


Best-fit 




4000 


0.46 


10 


28 


40 


1678 





Table 12 

Power Law Models: Best-Fit Summary 



Source 


P 






T 




Lint 


SiSRF 


Pm 






(xlO^CTO-^) 


(AU) 


(xloW) 




Le 






CB230 


1.9 


22 


100 


10 


0H5 


6.7 


0.3 


2.5 


IRAS 04166+2706 


1.6 


3.0 


100 


3.3 


0H5 


0.4 


0.3 


2.4 


IRAS 04169+2702 


1.5 


2.8 


100 


3.3 


0H5 


0.77 


0.3 


2.4 


IRAS 04239+2436 


1.8 


0.9 


100 


3.3 


0H5 


0.9 


0.3 


2.1 


IRAS 04248+2612 


0.8 


0.35 


100 


3.3 


0H5 


0.6 


0.3 


1.5 


IRAS 04295+2251 


2.3 


2.0 


100 


3.3 


0H5 


0.25 


0.3 




IRAS 04361-F2547 


1.8 


1.8 


100 


3.3 


0H5 


2.2 


0.3 


2.3 


IRAS 04381-F2540 


1.6 


1.8 


100 


3.3 


0H5 


0.53 


0.3 


1.7 


L1251B 


1.5 


12 


100 


7.2 


0H5 


10 


0.3 


2.3 



